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5^ . Abstract 

We present a detailed study of quantized noncompact, nonlinear S0(1,A^) 
^ . sigma-models in arbitrary space-time dimensions D > 2, with the focus on 

^ I issues of spontaneous symmetry breaking of boost and rotation elements of 

the symmetry group. The models are defined on a lattice both in terms of a 
transfer matrix and by an appropriately gauge-fixed Euclidean functional in- 
tegral. The main results in all dimensions > 2 are: (i) On a finite lattice the 
systems have infinitely many nonnormalizable ground states transforming 
irreducibly under a nontrivial representation of S0(1,A^); (ii) the S0(1, A^) 
symmetry is spontaneously broken. For D = 2 this shows that the systems 
evade the Mermin- Wagner theorem. In this case in addition: (iii) Ward iden- 
tities for the Noether currents are derived to verify numerically the absence 
of explicit symmetry breaking; (iv) numerical results are presented for the 
two-point functions of the spin field and the Noether current as well as a 
new order parameter; (v) in a large saddle-point analysis the dynamically 
generated squared mass is found to be negative and of order In V) in the 
volume, the 0-component of the spin field diverges as \/\nV, while S0(1, A^) 
invariant quantities remain finite. 



1. Introduction 



Noncompact nonlinear sigma-models occur in a variety of contexts. They are ubiqui- 
tous in the dimensional reduction of (super)-gravity theories, which provided the main 
incentive for the study of their quantum properties [1] - [7]. Motivated by structural 
similarities they were also used as a test-bed for renormalization and symmetry aspects 
of quantum gravity [9, 10]. The two-dimensional versions are in addition relevant for 
the theory of disordered systems and localization, see e.g. [12, 13, 14, 15, 16]. 

The most intriguing aspect of noncompact sigma-modcls is the apparent clash between 
symmetry and unitarity: the Lagrangian is invariant under a finite dimensional - hence 
nonunitary - representation of the group, while the physical Hilbert space (or at least 
a sizeable subspace of it) is expected to carry a unitary and hence infinite dimensional 
representation of the group, apparently not accounted for by the field content of the 
system. This is particularly puzzling in the vacuum sector, where in the 2-dimensional 
versions Coleman's theorem [17] seems to preclude spontaneous symmetry breaking even 
for a noncompact group. Indeed both perturbation theory and large N techniques 
typically expand around an invariant Fock vacuum in an indefinite metric state space 
[8, 1,5]. Its positive metric subspace, however, then carries no remnant of the original 
noncompact symmetry and looks more like that of a compact model. 

A recent detailed study of the 1-dimensional hyperbolic spin chain [18] showed how in 
that system the clash is avoided: there are infinitely many nonnormalizable ground states 
transforming under an irreducible representation of the group. On the one hand this 
entails that the symmetry is spontaneously broken at the level of (certain) correlation 
functions. On the other hand, by a change of scalar product to the one induced by 
the Osterwalder-Schrader reconstruction, the above representation rotating the ground 
states into each other can be made unitary. The price to pay is that the reconstructed 
Hilbert space is nonseparable and that the unitarity of the representation only extends 
to a 'large' but proper subspace of it. One of the goals of the present paper is to 
investigate the extent to which this picture of the 'ground state orbit' generalizes to the 
field theoretical case. 

More generally our focus is on issues of spontaneous symmetry breaking of non-compact 
(boost) and compact (rotation) symmetries. The starting point is a lattice construction 
of the models, using both the transfer matrix formalism and the Euclidean functional 
integral. In either case the infinite volume of the symmetry group requires modifications 
compared to the setting for a compact symmetry group. The transfer operator is no 
longer trace class even in finite volume and the functional integral needs to be gauge 
fixed. Two specific gauge-fixing schemes (a translationally invariant scheme in which 
the zero-momentum mode of the transverse spin fields is set to zero, and a fixed-spin 
gauge) are used, the first of which is convenient for numerical simulations while the 
second one allows one to relate the transfer matrix to the functional integral. Once 
properly defined (Section 2) the systems are studied by a combination of group theoreti- 
cal techniques (Section 3), numerical simulations (Section 4), and a large N saddle-point 
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analysis (Section 5). Our main results in generic dimensions D >2 are: 

• On a finite spatial lattice the noncompact models arc shown to have infinitely 
many non-normalizable ground states transforming irreducibly under S0(1, A^) - 
in sharp contrast to the unique ground state of the SO(l + iV) models. 

• Spontaneous symmetry breaking occurs in all dimensions D >2. 

As described, the symmetry breaking is surprising in dimension D = 2; a, case which we 
therefore investigated in more detail with the following results: 

• A new 'Tanh' order parameter is used to probe the spontaneous breaking of the 
boost symmetries, bypassing problems with the usual 'hysteresis criterion'. 

• Quadratic Ward identities for the Noether currents are derived (including finite 
volume corrections) and used to verify numerically the disappearance of explicit 
breaking of the boost and rotation symmetries with increasing volume. 

• Numerical results are presented for the two-point functions of the spin fields and 
of the Noether current, as well as for the 'Tanh' order parameter, which show 
spontaneous symmetry breaking. 

• In a large saddle-point analysis (starting from the gauge fixed functional inte- 
gral) the dynamically generated squared mass is found to be negative and of order 
l/{VlnV) in the volume V, the 0-component of the spin field diverges as y/\nV 
while S0(1,N) invariant combinations remain finite. 

In addition we point out certain subtleties, related once more to the noncompactness of 
the symmetry group, without attempting definite answers here. One of them concerns 
the inapplicability of standard theorems in D = 2 (Mermin- Wagner, and refinements 
thereof) to argue that the maximal compact SO(A^) subgroup singled out by the gauge 
fixing is not spontaneously broken; see Section 2.2 for a discussion. For any D > 2 
another subtle point is the reconstruction of a Hilbert space, a transfer operator, a nor- 
malizable ground state and a representation of the symmetry group commuting with it 
from the infinite volume correlation functions via an Osterwalder-Schrader reconstruc- 
tion; see Section 3.5. 

The paper is organized as follows: in the next section we introduce the ingredients of a 
lattice construction of the systems (transfer matrix and functional integral) and pose the 
questions we wish to address. In Section 3 we derive the a structural characterization 
of the ground states of the finite lattice systems in D > 2 and prove that whenever a 
thermodynamic limit exists it shows spontaneous symmetry breaking. Sections 4 and 
5 are devoted to the D — 2 model, and contain the Monte-Carlo study of the S0(1, 2) 
models and the large saddle-point analysis, respectively. Some technical material on 
the harmonic analysis of functions on the target space and the finite volume corrections 
to the Ward identities are relegated to Appendices A and B, respectively. 
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2. Lattice construction 



Wc consider the hyperbolic S0(1, N) nonhnear sigma-models with N > 2 defined on a 
/^-dimensional Euclidean lattice, A C Z^, with D = d + 1 > 2. The systems are defined 
on finite lattices with the thermodynamic limit A — > to be taken later on. We divide 
the lattice into time slices At = {{x,t) \ 1 < < Lg , // = 1, . . . , d} ~ {1, . . . , L^}'^ 
of |At| = lattice points. The Euclidean time t ranges from to — 1, so that 
|A| = LtLg is the total lattice volume. The dynamical variables ('spins') are denoted 
by n^;, a; G A; those in a given time slice are written alternatively as n^:, a; G A^ or as 
'iT-x,t: a:; G {1, . . . , L^. The spins take values in Hjv — {n E B}'^ \ n ■ n — +1, nP > 0}. 
The bihnear form (dot product) is a • & = a°6° — a^b^ — . . .a^b^ —: a^W — ab, with 
a — {a},. . a^). We take as our basic lattice action 

where fi is the unit vector in the positive /^-direction (with ii = 1, . . . ,D, and the 
boundary conditions specified later). Since n • n' > 1 for all n,n' G Hjv the action is 
normalized such that So[n\ > 0. 

We use the connected component of the identity SOo(l, N) (which preserves both sheets 
of the cone a ■ a = 0) thoughout. Slightly simplifying (and abusing) the notation we 
shall always write S0(1,A^) := SOo(l.A^) for it. The hyperboloid EItv can then also 
be viewed as a globally symmetric space S0(1, iV)/SO(iV) for any one of the maximal 
compact SO(A^) subgroups. We shall use the stabihzer group SO^(A^) of the vector 
— (1, . . . , 0) throughout. Concretely this amounts to a parametrization of the 
spins as n = (^, — Is), where > 1 is a noncompact variable and s G S'^~^ is 
a conventional compact spin. Note that this provides a global parametrization of Hjv. 
The invariant products entering the action then read 

Tlx ■ n,+A = Ux+f. - iil - l)'/'(e'+A - ■ 4+A , (2.2) 

We write 5*0 for the action in this parametrization. It can be viewed as that of 
a spherical S'^'^ sigma-model coupled in a non-polynomial way to the additional non- 
compact field The invariant measure dQ{ri) :— 2d^^^nS{n-n — l)9{'nP) factorizes 
according to 

//■oo r 
dn{n) = / dae - If^'-' / dSis) . (2.3) 
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2.1 Definition of the transfer matrix and functional integral 



The dynamics of the lattice system is defined in terms of the transfer operator T which 
transports lattice configurations from one time slice to the next. The square integrable 
wave functions ipln) = ip{nx, x G A^) depending on the spins in some time slice At 

form the Hilbert space L'^(H.^") with respect to the product of the canonical invariant 
measure. Since we usually keep Lg fixed we simply write for this Hilbert space once 
the number of spatial dimensions d is clear from the context. The transfer operator T 
acts on as an integral operator via 

(T^)(n) = fl[dn{n'^)rf,{n,n';l)^{n'), (2.4) 
xeAt 

Tp{n, n'; 1) = D^^^^ exp i -/3 [^^ ' < + \<^^ ' %+i + \^'x " <+i - 2 

The normalization constant -D^ jv is introduced for later use; it sets the overall scale in 
that < Tp(n,n'; 1) < Dp'^ for all configurations. The variables in (T'0)(n) can then 
naturally be associated with the time slice A^+i. Indeed upon iteration of (2.4) one 
obtains 

(TV)(n) = /"n dn{n'^)Tf,{n,n';t)ij{n') , (2.5) 

xeAo 

where we conventionally regard T* as a map from time slice Aq to time slice A^. In this 
interpretation the iterated kernel reads 



Tfs{n e At,n e Ao;t) = 0^,^^/ exp i | ^ ( ^ - ^ )n, ■ 

I H^D xeAo xeAt 
a;eAi,...,At_i ^ ij.,xeAo,...,At-i 



■x+fi 



(2.6) 



For t — Ld and periodic be in the x/j-direction the last expression clearly resembles the 
partition function for the action (2.1). The last integration over the variables Uxfi — 
nx,L2J however, would be divergent as the infinite volume of Mn gets overcounted. 

To see this more clearly note that the wave functions ijj{n) carry the following "diagonal 
action" of SO(l,Ar), 

p{A)iP{nx,xeAt)=ip{A-^nx,xeAt), AeSO{l,N). (2.7) 

In contrast to a lattice system with a compact symmetry group p invariant wave func- 
tions, i.e. those satisfying p{A)ip{ri) — ip{ri), for all A 6 S0(1, A^), do noi lie in the Hilbert 
space. This is because in the inner product on one of the integrations factorizes, the 
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infinite volume of Hjv gets overcounted and tlie norm of the wave function diverges. 
On an integral operator K with kernel K,{n,n') the group acts as X — > p{A)~^Kp{A) 
and thus as K.{n, n') — > K{An, An') on the kernels. In particular operators K whose 
kernels only depend on the invariants 

Tlx * '^y invariant. Importantly this holds for 

the iterated transfer operator, i.e. 

T*op = poT*, VieN. (2.8) 

Since the semigroup T*, t G N, describes the evolution of the system in Euclidean time 
Eq. (2.8) means that the dynamics is S0(1, N) invariant, as required. On the other hand 
it also implies that, although T is bounded, in contrast to the transfer operator of most 
other lattice systems in finite volume (as we shall see later) it is not trace class (see [18] 
for the 1-dim. case) . As a consequence correlation functions cannot be defined in terms 
of the usual expressions involving traces. The remedy is to ('gauge') fix the residual p 
symmetry by a variant of the familiar Faddeev-Popov procedure. To the best of our 
knowledge this gauge-fixing does not seem to have been taken into account in earlier 
studies, rendering the results somewhat formal. We now first describe this procedure 
and then outline the relation to the transfer operator. 

We used the following two gauge-fixing choices, both of which leave the stability group 
SO^(A^) of the vector = (1, 0, . . . , 0) intact: 

1. The noncompact global gauge freedom is eliminated with a translationally invari- 
ant gauge choice which sets the zero momentum mode of the transverse sigma 
fields fi :— a/^^ — Is to zero: 

^n, = 0. (2.9) 

x-eA 

In this gauge there is a nontrivial Faddecv-Popov determinant which comes out to 
be (^^ ^x)^! sec [19] in the compact case. The expectations of a general multilocal 
observable 0{{n}) then assume the form 

xexpj -,So[n] + Arin^n°}, (2.10) 

X 

where Zi{A, (3) is the partition function normalizing the averages, (1)a,/3,i = 1. 

2. Alternatively, the noncompact global gauge- freedom may be eliminated by freezing 
a single spin at an arbitrary site Xq to a conventional fixed value, typically = n' , 
where = (1, 0, . . . , 0) and G Aq. In this case there is no Faddeev-Popov factor 
and the expectation value of a general observable O is simply 

^^^^■^'^ ^ z;u^) I n ^({^» ' (2-11) 
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if we assume that the support of the observable does not include xq (otherwise 
an explicit 5{nxo,n^) factor has to be included). Although the choice of a specific 
site Xq would seem to destroy the translational invariance of the theory, the fact 
that this choice corresponds to a global gauge transformation implies that S0(1, A^) 
invariant observables are unaffected, and translational invariance still holds for such 
observables, provided of course this invariance is not explicitly broken by boundary 
conditions. We consider again periodic boundary conditions (be) and denote the 
resulting expectations by (C)a,/3,2- Invariant observables then should have the 
same expectations as with the translationally invariant gauge fixing, i.e. (0)a,/3,i = 
(C)a,/3,2, for O S0(1, N) invariant. 

We state without proof that the finite volume partition functions Zi and Z2 are well- 
defined, i.e. the gauge fixing is sufficient to render the integrals finite. Another inter- 
esting choice of be, in the case of the fixed spin gauge, would be periodic be in the 
spatial and free be in the temporal direction. In analogy to the 1-dimensional model 
the thermodynamic limit of these expectations should be expected to be different from 
each other, thereby revealing a peculiar kind of 'long-range order'. However to study 
the issue through numerical simulations would presumably require much larger lattices 
and a cluster algorithm. 

One can view ( )A,i3,i as hnear functional over the algebra of bounded observables C^, 
that is, continuous bounded functions O of finitely many spins with pointwise addition 
and multiplication and equipped with the supremum norm, O ^ As such they 

qualify as states in the statistical mechanics sense: < and for nonnegative O 

the expectation value is nonnegative. 

The fixed spin gauge with periodic be also allows one to make contact with the transfer 
matrix (2.6). For example 



/ 



II dnMS{n,„n^)Tp{n,n;Lt) ^ Z,{A,P) , (2.12) 



gives the partition function. Here Xq G Aq and 5{n, n^) is the invariant delta-distribution 
concentrated at n = with respect to the measure dQ{n). Similarly for the expecta- 
tion of a generic (noninvariant) observable (9(n^,%) located at x = {xi, . . . ,xd), y — 
(yi, . . .,1/0), one has 

{0{n^,ny))A^p^2 = 2 (\ g) [Yldfl{n^)6{n,=^^,n^) Y[ dfl{n^) Y[ dfl{ny) 

xTf^i^riz, n^] XD)0{n^, Uy) Tp{n^, Uy] yo - XD)Tp{ny, riz] Lt - yo) ■ (2.13) 

Because of the gauge fixing for finite Lt this is in general not invariant under time 
translations (nor, for that matter, under space translations). An important exception 
are S0(1,A^) invariant observables O, satisfying 0{Anx, Auy) = 0{nx,ny), for all A G 



6 



SO(l,iV). In this case Eq. (2.13) simplifies to 




(2.14) 



using the invariance of the integration measure and the convolution property for the 
kernels (2.6). Here translation invariance in the time direction (and trivially in the 
spatial direction) is manifest. Eqs. (2.13) and (2.14) generahze straightforwardly to 
observables O depending on more than two spins. 

Note also that any single SO(l,iV) transformation on an observable O can always be 
compensated by a change in the gauge fixing condition 



If be other than periodic ones were adopted the be would hkewise be "counter rotated" 
in (2.15). The issue of spontaneous symmetry breaking we wish to address in the fol- 
lowing, of course, asks for the invariance or noninvariance of the expectations under all 
of S0(1, A^) or a continuous subgroup thereof, with the gauge fixing and be held fixed 
and A ^ Z^. 



2.2 Ward identities — absence of explicit symmetry breaking 

In this context it is important to make sure that no explicit breaking of the symmetry 
is induced by the gauge fixing. In the thermodynamic limit one expects that the effect 
of a single fixed mode fades out, but experience with the ID model [18] shows that 
such expectations can be misleading. In addition, as simulations are done on a finite 
lattice a quantitative assessment would be useful. This can be done by deriving Ward 
identities expressing the S0(1,A^) invariance of all but the gauge fixing terms in the 
functional measure, with the latter giving rise to finite volume corrections of the 'naive' 
Ward identities. In this section we describe the principle of the derivation as well as 
the 'naive' form of the Ward identities which is dimension independent. In contrast the 
form of the finite volume corrections is dimension dependent and their determination is 
rather technical. For the case D = 2 (where the symmetry breaking is surprising) this 
is done in appendix B. Throughout this section the translation invariant gauge fixing 1 
with periodic be will be used, i.e. Eq. (2.10). 

In lattice models with a compact symmetry group the invariance of the functional mea- 
sure (including the Boltzmann factor) gives rise to Ward identities in a well-known way: 
implement a local symmetry transformation and expand the functional integral in pow- 
ers of the gauge parameter (s). Since the total response must vanish the coefficient of 



i = 1,2. 



(2.15) 
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each power must vanish, which gives rise to identities relating correlators of the Noether 
current to other correlators. In the case at hand the gauge fixing and the associated 
Faddeev-Popov determinant lead to a non-invariant overall measure. Nevertheless the 
impact of a local symmetry transformation can be computed and leads to modified Ward 
identities. 

We begin by fixing our conventions for the Noether current. It takes values in the Lie 
algebra so{l, N), and we normalize the components with respect to the previously used 
basis t"'^, < a < b < N, according to 

Jf{x) = pe'rix ■ = -Prix ■ = - «+^] ■ (2.16) 

Their two-point correlators can be decomposed into a transversal, a longitudinal, and a 
harmonic piece. This is conveniently done in Fourier space 

{Jfi^) Jf{y))A,p, =^ m E . (2-17) 

p 

where p runs over the dual lattice, = x'^m' % — 0; 1; ■ ■ ■ > ~ A* ~ 2- order 
not to clutter the notation we suppress the specifications (A, /3, i) on the right hand side 
(remember that i = 1,2 refers to the gauge-fixing adopted). The irreducible components 
Jx'ij)) (transversal), Jl'{p) (longitudinal), and Jnip) (harmonic) are picked out by 
acting with the corresponding projectors on J^^yip), see e.g. [20] for details. 

As mentioned earlier. Ward identities now arise from studying the response of a given 
expectation value under a local symmetry variation rix — > exp(Q;j;i"^)n-j:, x e A, per- 
formed on all spins. To get the response of the action we prepare (e^ := 77"" equals 1 for 
a = and —1 for a 7^ ): 

- ^(AA«.)'(e;,«+^ + 6Xn^+^) + 0{a^) ,(2.18) 

where A^a^. = a^+jx — and [{t"'^ynxY — —€bn'^6°^ — eaU^.S^'^ was used. For the change 
in the Boltzmann factor this gives 

exp{-/3^n^ -n^+^j — > exp{-/5 ^ • n^+^}| 1 + ^ A; J^''(a;) (2.19) 
+^ J2 «.a,A;jf(x)A;jf(y) + ^^(A^a.)2(e,<<+^ + e„<<^)+0(a3^ 

x,n;y,i' x,tx 

where A*^.Jjf'{x) = Jf'{x) — Jf'{x + jl). Using (2.19) the expansion of the Boltzmann 
factor to 0(q;^) is trivial. Since the product of the invariant measures ^[^'^^^^(^^ ~ ^) 
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is invariant even under local SO(l,iV) rotations, the only non-invariant terms in the 
functional measure come from the gauge fixing. Focusing on the invariant terms the 
total response under a local symmetry transformation must vanish. In principle the 
vanishing of the coefficients of each power in a gives rise to a new identity. 

For example, the 0{a) terms in the response of (npA,/3,i give rise to the following first 
order Ward identity 

(A; Jf npA,/3,i + 5,,y{{i''\Y)KAi + terms from gauge fixing = . (2.20) 

Replacing by with a generic (noninvariant) observable . . . , n,j.^^ a similar iden- 

tity arises where the correlator with A*^J^'' produces a sum of contact terms. We shall 
not pursue these first order Ward identities further: in Section 5 we shall verify in a 
large N analysis that (n")A,^,i diverges as |A| — > oo. One expects this to hold also at 
fixed A^. in which case already the example (2.20) shows that these first order Ward 
identities do not necessarily have an interesting thermodynamic limit. This very fact 
however is worth mentioning, because it shows how the confiict with Coleman's theorem 
[17] is avoided: the currents simply do not exist in the thermodynamic limit. 

More interesting is the second order Ward identity from the response of the partition 
function itself [20]. The vanishing of the 0{a'^) terms requires 

I J2 «.a,(A;jf A;jf)A,^, + |^(A^a)2(6,£;'^ + 6„£;'') 

x,iJ,,y,i/ x,fi 

+ terms from gauge fixing = , (2-21) 

where E"' := -E^/Ji {''^x''^x+fi)^,l3,i action link variable. The terms induced by 

the gauge fixing can in principle be computed exactly. They are expected to die out as 
A lP ^ but the precise form of the correction terms is cumbersome to compute. 

As the case D = 2 is of particular interest the derivation of the finite volume corrections 
is detailed in appendix B. The extra terms induced by the gauge fixing then turn out 
to be of order 0(ln|A|/|A|) in the hmit of large volumes |A|. Converting (2.21) into 
Fourier space the longitudinal part J'{^{p) of the current two-point function appears. 
The resulting Ward identity generalizes that in the compact models [20] and reads 

Jl\v) = -ii{e^E^ + e^E'>)-r0{^^^, Vp , a<6, (2.22) 

On account of the invariance of the vacua under the maximal compact subgroup singled 
out by the gauge fixing (see Section 2.3) one expects that E'^ > — . . . — E^ > 0, so 
that only two distinct cases arise: 

Jl'^{p) = +2f3E^ + o(^-^^ , rotations, 

Jl\p) = (3{E''-E') + 0^-^y boosts. (2.23) 
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All quantities in (2.22), (2.23) of course depend on the specifications {A,(3,i). The 
inequahty > a 7^ 0, follows from nP > \n\. Combined with the trivial identity 
TT-x ■ iT-x+p, > 1 one gets the stronger bound — NE^ > 1. 

The individual E^ cannot be expected to have a finite thermodynamic limit. In Section 
5.1 we verify that in the large N expansion both E^ and E'^ diverge logarithmically 
with the volume, according to ~ ^ In |A| and E^ ~ -i:^ In |A|, where A = N/ jS. In 
contrast the invariant combination E^ — NE^ approaches the finite constant 1 + A/4. 
For the Ward identities (2.22) therefore only the invariant combination is assured to 
have a finite thermodynamic limit. Nevertheless the Ward identities for the individual 
components arc useful to test quantitatively the degree to which the boost/rotation 
symmetry is restored on a finite lattice (as far as the dynamics is concerned). Since the 
current correlator (2.17) and the action link variables E"- are independently measurable 
quantities in a Monte-Carlo simulation, validity of the identities (2.23) also provides a 
good test on the simulations for given lattice size and boundary conditions. We report 
the results of such a test in Section 4.1. 



2.3 Spontaneous symmetry breaking and 'Tanh' order parameter 

Even the very notion of spontaneous symmetry breaking in the noncompact models 
requires a little thought. The conventional analysis of spontaneous symmetry breaking 
asks if there is a local observable having a noninvariant expectation value if we either 

(a) fix symmetry breaking boundary conditions and then take the thermodynamic 
limit, or 

(b) add a symmetry breaking term like a magnetic field h to the action, take the 
thermodynamic limit and then turn the symmetry breaking term off. 

In the second picture spontaneous symmetry breaking amounts to a hysteresis effect. In 
a model with a compact symmetry group then the one-sided derivatives /i 0^ and 
/i — > 0^ exist, but are different. This way of looking at spontaneous symmetry breaking, 
however, does not readily generalize to the boost symmetries in the noncompact sigma- 
models because the field h has to serve double duty - as a regulator and as a probe for 
symmetry breaking. For invariant observables it is clear that a nonzero field is needed 
in order to (potentially) produce a normalizable measure even in finite vohimc. For 
noninvariant observables coupling to a magnetic field may or may not not render the 
finite volume expections finite. Indeed, a typical coupling would add a term of the form 
to the action (2.1). However for /i < then already the finite volume averages 
fail to exist. The /i — > 0+ derivative is expected to be convergent for D > 3 and divergent 
for D = 2. 

Indeed, since the first version of this paper was posted, an interesting result by Spencer 

and Zirnbaucr appeared [21], in which it was shown that in D > 3 the expectations 
(n°)A,/3,/i defined without gauge fixing and with a positive magnetic field h can be 
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bounded by a constant (independent of h and |A|) for all (3 > 3/2 and \A\h > 1. 
Thus a 'one-sided hysteresis criterion' here signals spontaneous symmetry breaking in 
the thermodynamic limit. 

The case D — 2 will be discussed in more detail in section 2.4; in this case even the 'one- 
sided hysteresis criterion' is expected to fail. In D = 2 many authors found that 
diverges in the thermodynamic limit, based on an (un-gauge fixed) large expansion. 
This amounts to some vestige of the large fluctuations that are responsable for the sym- 
metry restoration in compact and abehan models. However since {A'nP),A e 80(1,^^), 
diverges likewise one can only conclude that the symmetry breaking cannot be seen on 
this particular observable. 

The approach adopted here is somewhat different. The gauge fixed functional integrals 
(2.10) and (2.11) provide a complete definition of the systems in finite volume, both for 
invariant and for non-invariant observables. The regulator (gauge fixing) is decoupled 
from whatever probe is used for the symmetry breaking. Spontaneous symmetry break- 
ing can then be discussed without appeal to a 'one-sided hysteresis criterion' and for all 
D > 1. The criterion we propose is simply that there exist noninvariant observables O 
for which the thermodynamic limit exists and for which 

lim {0{An))A,p,i ^ lim (0(n))A,^,i , for some A e S0(1, TV) . (2.24) 

Wc should remark that (2.24) for a boost A signals also breaking of the compact subgroup 
obtained by conjugating SO^A^) with A. Spontaneous symmetry breaking then basically 
follows from the nonamenability of the group S0(1, A^). For convenience we recall the 
definition here: 

A Lie group G is called amenable if there exists a left invariant positive linear functional 
("a mean") on Cb{G), the space (and commutative C*-algebra with unit) of bounded 
continuous functions on G equipped with the sup-norm. All abelian and all compact Lie 

groups are amenable. Conversely, G is called nonamenahle if no such mean exists. All 
noncompact scmisimplc nonabelian Lie groups arc known to be nonamenable. 

In the present context the nonamenability of S0(1,A^) implies that there has to be 
bounded continuous functions of one spin, say at the origin, whose infinite volume ex- 
pectation values are not invariant under the group. The precise form of this result is 
described in Theorem 3 of section 3. 

In [18] we identified for the ID model a useful example of such a function, the so-called 
'Tanh' order parameter. As explained in Section 2.2 the 'one-sided hysteresis criterion' 
to describe spontaneous symmetry breaking cannot readily be used for D < 2. The 
'Tanh' order parameter, on the other hand, docs not require the introduction of an 
external field; the gauge fixing or the boundary conditions single out the direction of 
symmetry breaking and the maximally compact subgroup SO^(A^) that remains unbroken 
is the stabihty group of . This construction readily generalizes to all > 1: 
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For a spacelike unit vector e we define 

Te(n) := tanh(e ■ n) , e ■ e = —1 , 

:= / dpL{A)T,{An) (2.25) 

JS0\N) 

= / dfi{A) tanh U^/^^ - gVC'-leo • As) , 

JsO\N) ^ ^ 

where dfj,{A) is the normahzed Haar measure on SO^(A'"). After the group averaging the 
observable only depends on ,^ := n^^ ■ n and q := \/n^ ■ e + 1. Here we parameterized n 
and e as n = a/^^ — Is), s^ = 1 and e = (^g^ — 1, geo), = 1- This observable 
of course remains finite for A — >■ Z-^ even if {pP)K,i3,i diverges. More importantly it is 
designed to be a good indicator for 'spontaneous' symmetry breaking already in finite 
volume. The criterion (2.24) for spontaneous symmetry breaking becomes for all D >1: 
(Te(n))oo„fl,i 7^ {Te{An))oo,i3,i for some A G S0(1, A^). Since by Section 2.3 the finite vol- 
ume average in itself effects the S0\A^) average this is equivalent to T[q) := {Tq{n^))oo,f},i 
having a nontrivial dependence on q. Clearly \T[q)\ < 1 and T(l) = 0, by the SO^(A^) 
invariance. 

Typically a nonzero value for {Tf,{n))A,i3,i at some g > 1 is numerically easy to detect. 
In order to view this as a signal for spontaneous symmetry breaking one has to exclude 
that this value decays to zero as A — > Z^. Since by a convexity argument one expects 

{Te{n'))A,p,i > T,((n°)A,^,,) > T,(sup(n°)A,^,) , (2.26) 

every nontrivial {n^)A,i3,i will thus provide a lower bound on the measured (Te(n))A,/3,i > 
0, which therefore cannot decay to zero as A — > Z^. For D = 2 we shall find later in the 
large N limit that T{q) is in fact a strictly increasing function of q approaching 1 for 
g — > oo. Specifically one has Tq(^) ~ tanh(ri-\/g^ — 1), with n given by Eq. (5.6) below. 
For N = 2 the same monotone increasing behavior is found in numerical simulations, 
see Section 4.3. 



2.4 Unbroken S(y{N) invariance in £> = 2? 

In two dimensions an additional subtlety arises from the Mermin- Wagner theorem [22] 
and its refinements [23, 24]. Whether in the fixed spin gauge or in the translation 
invariant gauge, the system has a residual SO^(A^) invariance and can be viewed as a 
O(A^) vector model with fluctuating length of the spin vectors. In D = 2, at flrst sight 
it may seem obvious that this compact symmetry cannot be spontaneously broken, due 
to the mentioned theorems. 

On closer inspection, however, the situation is not quite as simple. The above mentioned 
theorems on the absence of spontaneous symmetry breaking cannot really be applied. 
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because some technical conditions for their apphcabihty are not fulfilled. The first one 
is the condition that the second derivatives of the interaction with respect to the group 
parameters have to be uniformly bounded over the configuration space of the spins, 
which fails as a consequence of the noncompact nature of the latter; c.f. (2.2). The 
second condition is that in the thermodynamic limit we have to have a Gibbs measure 
on the configuration space, whereas in fact our infinite volume state is not a measure, 
but only a more general mean (for more detailed discussion of this see [18]). 

The symmetry breaking be in option (a) of section 2.3 must break the noncompact 
symmetries and hence amount to something similar to the fixed spin gauge. Also here 
some potential pitfalls arise, which wc illustrate now for the S0(l,2) model. If one 
looks at individual configurations of the S0(1, 2) model in the fixed spin gauge at weak 
coupling (specifically, /5 = 10 and nxo=o = 0, say), one finds that contrary to the 
naive expectation expressed above, all spin vectors fix seem to point roughly along the 
same direction: the system appears to have acquired a large spontaneous magnetization! 
Writing momentarily ( • )i for ( • )a,/3,i with Lg — Lt —: L, we find that the 'average 
magnetization' defined as 



does not vanish, rather seems to increase with L: On a 32^ lattice it is 0.7197, on 64^ 
it is 0.7277, and on 128^, 0.7361 . The reason for this behavior becomes clearer if one 
considers the 0-component of the spin: while it is fixed to unity at the origin, it grows 
logarithmically with the distance. But large zero components necessarily imply large 
"spatial" components {as n ■ n — 1), and consequently a large ferromagnetic coupling 
between neighboring spins. 

To understand this phenomenon in more detail, it is useful to look at the zero-curvature 
limit of the S0(1, 2) model, which is just the Gaussian model of a two-component mass- 
less free field, defined by the lattice action 



In this model the fixed spin condition means that the spin at the origin no is fixed to 0. 
Furthermore, we can compute the counterpart of (2.27) analytically. First one has 




(2.27) 




(2.28) 




k 



(2.29) 



where the well-known function D is given by 



1 



exp 



2iTil-x 




(2.30) 
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with the Yl' running over /i, /2 = 0, 1, . . . L — 1 but h = I2 = omitted. Using this, we 
obtain for the square of the numerator of (2.27) 



^((E-^)'), = ^E^(-)' (2-31) 



X X 



and for the square of the denominator 

X X 

It is apparent that in this model one gets M = l/\/2 = .7071, independent of the lattice 
size L. The closeness of this number to the numbers quoted above is striking. The small 
difference between the numbers is due to the curvature of the target space of the S0(1, 2) 
model, which becomes relevant as the spins fluctuate further from the fixed spin at the 
origin. This explains why the difference grows with growing L. On the other hand, by 
increasing (3 the curvature should become less important; we checked this by measuring 
M at /3 = 40 on a 8^ lattice and found M — .707 in agreement with these expectations. 

We conclude that the apparent magnetic ordering of the lattice is due to the fact that the 
spins fluctuate very far from the origin where no — 0, and these excursions necessarily 
take place in a certain direction. If we think of the spins of the noncompact model as the 
on-mass-shell momentum vectors of a unit mass particle, these vectors are constrained 
by the action to be such that neighboring particles on the lattice are roughly collinear at 
weak coupling. In the fixed spin gauge, the fixed vector at the origin then corresponds 
to a particle at rest, surrounded by "nonrelativistic" neighbors. Far from the origin, 
the particles become highly relativistic, collinear, and with a local center of momentum 
frame highly boosted relative to the rest frame of the spin at the origin. This global drift 
of the center of momentum frame can be prevented by a choice of gauge: indeed, this 
is precisely the role of the gauge-fixing condition in the translationally invariant gauge, 
where the spins are described in a frame in which the total spatial momentum vanishes. 
The fixing of only a single spin is insufficient to arrest the gradual drift of the global 
center of mass frame in the large volume limit: this phenomenon happens likewise in our 
model and in the two-component massless free field. Nevertheless one would not ascribe 
spontaneous symmetry breaking to a Gaussian model. 

In fact, this ordering of the lattice is really a gauge artifact: Obviously, in the translation 
invariant gauge J2x^x vanishes and so does the magnetization (2.27). We have checked 
that by boosting the configurations from the fixed spin gauge to the translation invariant 
gauge, the magnetization disappears - instead of one dominant direction of the spins we 
find domains of different spin orientations. 

The upshot is that in the discussion of spontaneous symmetry breaking viewpoint (a) 
should be adopted. The symmetry breaking be must break the noncompact symmetries 
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and hence amount to something similar to a gauge fixing of a single spin leaving a 
maximal compact subgroup, here SO^(iV), intact. In this gauge we could not find any 
local observable that has a thermodynamic hmit and shows breaking of the SO^(A'") 
symmetry, and by analogy with the Gaussian model discussed above, we do not think 
that such an observable exists. The translation invariant gauge fixing is sometimes 
more convenient but should lead to the same conclusion. In summary, although the 
usual theorems do not apply, we expect that for local observables in the two-dimensional 
model 

{p{A)0)^,fS,i^{0)ooAi^ i = l,2, yAeSO\N). (2.33) 

Here oo refers to a two-sided thermodynamic limit where A — > Z^. We shall offer some 
further comments on (2.33) in the conclusions. 



15 



3. The ground state sector 

Spontaneous symmetry breaking, as described above on the level of correlation functions, 
only precludes the existence of an invariant ground state, without saying much about 
the set of possible ground states and the possible action of the symmetry group on 
them. In this section wc present a hamiltonian analysis of the ground state sector of the 
lattice systems in a finite spatial volume, both for discrete and for continuous Euclidean 
time, that results in a very concrete description of the 'ground state orbit' mentioned 
in the introduction. The discussion, though limited to systems of finite spatial extent, 
is valid for S0(1, N) sigma models in arbitrary (spatial) dimensions. An outlook on the 
thermodynamic limit via the Osterwalder-Schrader reconstruction is given in Section 
3.5. 



3.1 Transfer operator in the Schrodinger representation 

To address structural issues the transfer operator in the Schrodinger representation is 
useful. We begin its construction by describing the infinitesimal form of the S0(1, A^) 
representation p in (2.7). Recall that we write S0(1,A^) for SOo(l,A^). Then A e 
S0(1, N) if and only if A preserves the bilinear form a ■ b = a^lP — a^b^ — ... — a^b^ 
on R^'^ and both sheets of the cone a ■ a = 0, and has unit determinant. In matrix 
components the first condition becomes A^'^rjcdAf^'^ = rjab, where r] = diag(l, — 1, . . . , — 1), 
while the second condition amounts to > 0. For elements t of the Lie algebra 
so{l,N) the defining relation reads ta ■ b + a ■ tb — 0, or t^^'^ricb + Vadtb^ — 0. An exphcit 
basis is {t"'^)^'^ = — S^rj"-^, < a < b < N. Consider the 1-parameter subgroups 
R e s I— > exp{st'^^}, generated by these basis elements and set 

Pif-'mn) = x e A,) [^^ = ^^11 ' 

* * xeAt 

with p^(t«'') = -(rV)^^ , xeAf (3.1) 

The differential operators Px{t"-^), < a < b < N, generate commuting copies of the Lie 
algebra so(l, A^) at each site: 

[p,{e') , pyit"')] = d^y [r^"V.(t'') - v'^'pAt''^) + v'^Pxin - v'^'Pxin] . (3.2) 

The normalization of p{t"'^) (generating identical infinitesimal transformations in all 

variables) is adjusted such that they satisfy the same algebra as the t°^. The quadratic 
Casimir of the local algebras coincides at each site with minus the Laplace-Beltrami 
operator on Hjv and is given by 

C, = -AH-= Yl PAf'^aa'Vbb'Pxit"'"'). (3.3) 
a<b,a',b' 
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We shall mainly need the following property of A'^'^ (omitting the site index mo- 
mentarily): the spectrum of — A''^^ is purely continuous and is given by the interval 
^{N — iy + u!^, cu > 0. Several complete orthonormal systems of improper eigenfunc- 
tions arc known, see Appendix A and e.g. [25, 26, 27]. In representation theoretical terms 
this corresponds to a decomposition of the quasi-regular representation on L^(HlAr) into 
a direct integral of unitary irreducible representations tt^^ of the principal series (where 
unitarity of the representation refers to an induced inner product; see e.g. [25], Vol.2, 
Section 10.1.4). For the representation spaces one has 

/■e 

L^(HAr) = / dou i^Nioj) Cn{oj) , (3.4) 

with an absolutely continuous spectral measure iin{ijj) duj given in (A. 3) and each irre- 
ducible representation occurring with unit multiplicity. Note that the (unitary) singlet 
representation is not contained in the decomposition (3.4). 

Consider now the integral operator T with kernel 

t(i{n-n';l) = D~^^exp{p{l - n ■ n')} , 



O iV-1 

where K^{z) is a modified Bessel function. The kernel of the iterated operator T^, a; e N, 
is denoted by t^(n • n'; x). The normalization is such that 



/ 



dQ{n)tp{n-n']x) = 1 , Vn' e H^v , a; G N . (3.6) 



T can be shown to be a bounded selfadjoint operator on L'^iM^) (see [18] for N = 2). 
The spectrum of T is absolutely continuous and can be computed exactly (see Appendix 
A) . The spectral values come out as 



= j^^'^fi , a;>0. (3.7) 



They are smooth even functions of u with a unique maximum at o; = 0. Although 
real and strictly bounded above by ii'o(/3)/-ft'(Ar-i)/2(/3) < 1 they are positive only for 
< u < u;+{P), where u;+{P) increases with /3 hke co+{P) ~ /3-|-const P^^^. For u; > u;+{/3) 
the behavior is oscillatory with exponentially decaying amplitude. Positivity however is 
restored in the (naive) continuum limit for the Euclidean time: introducing momentarily 
the lattice spacing a, continuum times r = xa, as well as a coupling = l/(/3a), one 
finds 

hmJA 1 .(a;)]^ =exp|-r^f^^^^ + a;2)l . (3.8) 



a— >0 g'-^a ' I 2 
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This allows one to make contact with the heat kernel, i.e. with the integral kernel 
ex'p{^A^^){n,n') of the exponentiated Laplace-Beltrami operator. Prom (3.8) one 
expects 

2 

exp ( ^A^A in. n') = lim t i fn • n'; [-1 V (3.9) 
\ 2 / a^o ^ V LaJ / 

where [x\ denotes the integer part of a; G M. On the one hand this can be shown to lead 
to the correct path integral (see Eq. (11.30) and Appendix D of [26]). On the other hand 
one can insert the spectral resolution for T to obtain that of the heat kernel. Then the 
spectral values X^^n of T are simply replaced with their continuum counterparts (3.8), 
see (A. 13), (A. 15). Massaging this integral further one can show the equivalence to the 
usual expression for the heat kernel on Hjv, quoted in Appendix A. For our purposes a 
crucial property is the strict positivity 

2 

exp (^^A'"^)(n,n') > 0, Vn,n'eHjv. (3.10) 

For finite lattice spacing we can define T in the Schrodinger representation through its 
spectral resolution. That is, 

T = Vjv(-A=-), with 
W«) ■■= \p,N{^s-{N-lf/A) , s>h,N-lf. (3.11) 

These results for T directly carry over to the 'kinetic' part of the transfer operator Tq 
which we define through its integral kernel 



'P. 

xeAt xeAt 



with the same normalization constant as in (3.5). In the Schrodinger representation this 
gives 

To = n Viv(-A=-) . (3.13) 

xeAt 

In particular Tq has absolutely continuous spectrum given by 



<7(To) = (T,.e.(To) = { n Vjv(c^x) I > o} , (3.14) 



xeAo 



where the overbar refers to the closure in R. Its supremum, i.e. the spectral radius g{To) 
of To equals q{Tq) = jv(O)'^^. On the other hand To is clearly symmetric with respect 
to the inner product and it is also bounded in the norm, since T is. It follows that 
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To extends to a unique selfadjoint operator on L^. As such its norm coincides with 
its spectral radius: ||To|| = ^'(To) = A/3^Ar(0)^^. The improper eigenfunctions of Tq are 
of course just direct products of those of T (see Appendix A) and thus are manifestly 
non-normalizable with respect to the norm. 

To proceed we introduce (unbounded) multiplication operators x e {1, . . . , Lg}'', 
such that for any function V{n) on H^y'' (in the Schwartz space of rapidly decreasing 
smooth functions, say) V{n)%l){nx.i x e A^) = Vinx, x e A.t)^l){nx1 x e A^). Defining the 
'potential operator' 



one readily verifies from (3.12) and (2.4) the following expression for the transfer operator 
in the Schrodinger representation: 



It is again bounded and symmetric with respect to the l? inner product and hence 
defines a unique selfadjoint operator. The norm satisfies 



This means T is a contraction satisfying (V',T*'0) < (V','0), for all V e ^^ ^ e N. 
In particular limt^oo('^, T*-?/') = 0. The norm ||T|| is of particular interest because for 
bounded selfadjoint operator like T it coincides with the spectral radius and hence can 
be thought of as "the exponential of minus the ground state energy". In section 3.2 and 
3.3 we shall try to narrow in on the associated (improper) cigcnspace. 

Of course in contrast to Tq neither the spectrum of T nor its eigenfunctions are analyt- 
ically accessible (except, perhaps, for very small L^). However there are some generic 
properties which the transfer operator of most lattice systems in finite volume has, but 
which one should not expect T to have, taking the known properties of Tq as a guide- 
line. Since one is very much tempted to tacitly assume them, we list the expected 
'non-properties' here: 

• The eigenfunctions of T should not be expected to be normalizable with respect 
to the l? norm. 

• T should not be expected to be a positive operator; for discrete Euclidean times 
then no lattice Hamiltonian, — In T, exists. 

• The eigenspaces of T (considered for instance as spaces of smooth bounded func- 
tions) carry representations of S0(1, A^) (the restrictions of p) which by the first 
point cannot be expected to be unitary with respect to the norm. 




(3.15) 



T = exp { - ^-V{n)] To exp { - ^y(ri)} . 



(3.16) 




(3.17) 
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To the first point one can add that if there was a normalizable ground state it would 
necessarily have to be noninvariant, by the remark made after Eq. (2.7). The second 
point is purely technical and could be avoided by working with instead of T. Alter- 
natively one could start with a different lattice action: the heat kernel action would be 
an obvious choice. For our purposes the most natural way to ensure the existence of a 
Hamiltonian is to take the continuum limit in Euclidean time at the outset. This limit 
exists as we shall argue now. 

To this end we assign a coupling /?2 to the kinetic part in T and a coupling /5i to the 
potential part and write T/j^ for the result. The k-th. power of this operator is given 

by 

TJ... = e-^^(To,,,e-/^^^)'e^^ (3.18) 

in accordance with the integral kernel (2.6). We introduce the lattice spacing a in the 
Euclidean time direction, writing t = ka for the continuum time, and set P2 — 
Pi = For large P2 one has T ~ exp(2^A''^^), so if we heuristically replace T by this 
heat kernel, we are led to consider instead of (3.18) the sequence 

e te^i^x^x g \ g3fc^ ^ (3.19) 
which is recognized as Trotter approximant (for k — > oo) of 



exp 



where the operator in the exponent is interpreted as the self-adjoint operator given by the 
form sum of — y X^x ^-r ^ both operators are unbounded but positive, 

actually Kato's strong Trotter product formula [28] and its refinements [29] could be 
applied to show that (3.19) converges strongly and even in trace norm to (3.20). But 
since we used the heat kernel approximation to Tq above in an informal way, this does not 
lead to a rigorous proof of the Hamiltonian limit. Probably with more work this could be 
done, but we do not really need this here; we simply take the semigroup R+ 9 r i— > T^, 
where 



-rH 

,2 



H = J]A^ + ll-, (3.21) 

xeAt ^ 

as the definition of the continuum dynamics. The essential self-adjointness of H on the 
space of smooth functions of compact support can also be seen directly, using the results 
of [30, 31]. Tt is a strongly continuous contraction semigroup. Both the semigroup 
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T-r, T > 0, and its generator HI commute with p. From (3.21) and (3.17) one infers the 
bound ||T^|| < exp(-rLf c/^^jV- 1)2/8), so that 

a(H) > — (iV-l)2 . (3.22) 
8 

Since EI is an unbounded but manifestly positive operator one could now search for 
a ground state in the usual way. Technically it is more convenient to work with the 
bounded positive operator or the semigroup T^, r > 0. 



3.2 Existence of positive ground state wave functions 

We begin by showing that and have no normalizable ground states. Here the 
concept of a positivity preserving or positivity improving operator T is useful [32] . For 
convenience we recall the definitions. A nonzero function & L"^ is called positive if 

ip{n) > almost everywhere (a.e.) (that is, outside a set of measure zero in with 
respect to the product of the invariant measure on Hat) and strictly positive if ip{n) > 
a.e. Then T is called positivity preserving if (T^)(n) > 0, a.e. and positivity improving 
if {Tip){n) > 0, a.e. The latter is equivalent to (■0i,T''02) > for all positive ipi,ip2- 

The classic use of the positivity improving property is to establish the uniqueness of a 
ground state once it is known to be normalizable, see [33] or [34] for an application to 
gauge theories. Here the argument works somewhat differently: the positivity improving 
property entails that a groimd state cannot be normalizable with respect to the norm. 

Wc first show that both and the semigroup T,- are indeed positivity improving. For 
T*, t G N, this is obvious because of the strict positivity of the kernels 7g(n, n'; t) in (2.4) 
and (2.6), rendering all matrix elements with strictly positive functions strictly positive. 
To show that T,- is positivity improving we use the the path integral representation of 
this semigroup which is possible because of the strict positivity of the heat kernel (3.10,). 
From [26, 35, 36] one infers 

Tr{n,n')^ J ex.p(^- V{uj{t)d?jdiio{uj) , (3.23) 

where djiQiuj) is the measure describing Brownian motion on with starting configu- 
ration n and end configuration n' and V is (up to a trivial normalization) the potential 
defined in (3.15). The paths uj are continuous a.e. with respect to the measure d/ioi^^), 
and for this reason the integrand is strictly positive a.e., so strict positivity of the kernel 
T^(n,r7/) follows, i.e. is positivity improving. Clearly the T^, r > 0, are also posi- 
tive operators while for the transfer operator it convenient to work with the manifestly 
positive square T^. 
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To proceed we recall from (see [32], Section XIII. 12) the following general result: if 
||T|| is an eigenvalue with a (proper) eigenvector the eigenvector is nondegenerate 
and can be chosen strictly positive. The latter statement corresponds to the folklore 
that a (normalizable) ground state wave function does not have 'nodes'. Applied to the 
case at hand an immediate consequence of this result is that neither nor T,- have 
a normalizable ground state. This is because such a ground state would have to be 
unique, and thus (since and commute with p) would have to be a singlet under 
the S0(1, A^) action p. However we already know that p invariant wave functions are 
never normalizable because of the 'overcounting' of the group volume, and thus arrive 
at a contradiction. In other words, for the noncompact sigma-models ||T^|| ^ o'pp(T^), 
||T^|| ^ crpp(T^), where a^p denotes the pure point spectrum. We are thus faced with 
the unusual situation that ||T^|| and ||Tt-|| must lie in the continuous and hence in the 
essential spectrum of the corresponding operators. In fact [37], and have only 
essential spectrum 

a{T^) = (Jess(T') , a{Tr) = a^{Tr) ■ (3.24) 

Recall ([38], Section VII. 3) that for a bounded selfadjoint operator T the spectrum 
cr{T) decomposes into two disjoint sets, the discrete spectrum adisc{T) and the essential 
spectrum aess{T). where acss{T) is a closed subset of M. In terms of the spectral projectors 
Pj this amounts to the distinction: A G crdisc{T) iff the range RanP/ of Pj is finite 
dimensional for some open interval / containing A, and A e o'ess{T) otherwise. Weyl's 
criterion states that A e cr(T) if and only if there is a family of normalized vectors {ipn)neN 
such that hm„_,.oo 11(7" — A)'0„|| = 0. Further A e cress{T) if and only if the vectors ■0n 
can be chosen orthogonal, so that their weak limit vanishes, i.e. lim„_»oo('0, ■^n) = for 
all jjj e L\ 

Note that does also not have ground states in the weak sense, i.e. vectors i/Joo ^ L'^ 
satisfying {ijj, (T^ — ||T^||)-(/'oo) = 0, for aU ip G L'^. This is because such a ipoo would also 
be a normalizable ground state in the ordinary sense. All this of course applies to T,- as 
well. 

Definitions: In the following we denote by T a transfer operator without strong or 
weak ground states. By a transfer operator we shall mean a bounded selfadjoint operator 
that is positive as well as positivity improving (and possibly subject to some subsidary 
technical conditions). In this situation one will naturally search for weak * ground states 
of T, i.e. solutions of (0, (T — ||T||)f2) = for all in a suitable function space, where fl 
is a vector in the dual space. Specifically we take (f) & and call fl G L°° a generalized 
ground state of T if (0, (T — ||T||)Q) = for all cf) G L^. The set of generahzed ground 
states forms a linear subspace of L°° which we call the ground state sector Q{T) of 
T. The existence of generalized ground states which are moreover strictly positive L°° 
functions is guaranteed by a general result [37] , a special case of which we describe here. 

Let be a locally compact space and Q a regular cr-finite Borel measure on it; let 
M. X M. 3 (m, m') h- >• T(m, m') — > M+ be a function that is symmetric, continuous and 
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strictly positive, i.e. T(m, m') > a.e. We also assume 

sup / dfl{m') T (m, m') < oo . (3.25) 

m J 

The last condition is sufficient (but by no means necessary) to ensure that T defines 
a bounded operator from to for 1 < p < oo; see [39] p. 173 ff. The operator 
norm ||T||lp^lp = sup||^||^^j^||T0||j, is bounded by the integral in (3.25) and coincides 
with it for p = l,oo. Positivity of the kernel entails that T is positivity improving. 
Positivity of the operator (that is, of its spectrum) is not automatic. However if it is 
not satisfied we can switch to and the associated integral kernel, where positivity is 
manifest. Without much loss of generality we assume therefore the kernel to be such 
that T is positive. As a bounded symmetric operator on the integral operator defined 
by T(m, m') has a unique selfadjoint extension which we denote by the same symbol T. 
The kernel of will be denoted by T(m, m'; L) for L e N. In this situation T and all 
its powers are transfer operators in the sense of the previous definition. 

To prove the existence of generalized ground states as defined above, we need a technical 
assumption, namely that there exist an m* £ M (an extremizing configuration) and a 
subsequence {Lj)j^^ C N such that 

r(m, m; Lj) <T{m^,m^; Lj) , Vm e Al and Vj e N . (3.26) 
The existence of generahzed ground states is then guaranteed by 

Theorem 1: Let T be a positive integral operator with kernel T(m, m') satisfying the 
conditions hsted above. Let E M. he as in (3.26) and set 

%(vmj := . ^ rJr* — r , Vm,{m) := T{m,m^) , (3.27) 

where ( , ) is the inner product on L^. Using the assumption (3.26) it is shown in [37] 
that the L°° norms of flj{vm^) are bounded uniformly in j. Therefore there exists a 
subsequence {jk)keN such that the weak * limit 

^m, w*- lim Qj^{vm,) , (3.28) 

fe^oo 

exists; because {vm,,Qj{vm,)) = 1 this limit docs not vanish. It is a strictly positive 
function in and a generalized ground state for T, i.e. 

(0, (T - ||T||)Q„J = forall0eL^ (3.29) 

Though our present proof requires (3.26) we expect that the conclusion of the Theorem 
remains valid for a much larger class of transfer operators, which arc not necessarily 
integral operators, and where in particular the condition (3.26) can be dropped. 
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In the case at hand, all but property (3.26) are manifest for our transfer operator T^. We 
conjecture that satisfies (3.26), in fact with an extremizing configuration e H^" 
where all spins equal . In section 3.4 we present in Theorem 2 a stronger result for 
which in particular entails the existence of generalized ground states in the above sense. 

Either way, the transfer operator of the S0(1, A^) nonlinear sigma model possesses 
strictly positive generahzed ground states. Based on Theorem 1 (and the conjecture) 
it comes parameterized by a preferred configuration n* G H^", i.e. G L°°{M,]^). 
This then gives rise to an entire orbit {VLau^ i ^ ^ S0(1, A^)} of strictly positive gener- 
alized ground states. In compact models the counterpart of this orbit would be trivial, 
i.e. would consist simply of the one-dimensional unitary representation ilAn* — > for 
all A G SO(l-l-A^). It is a remarkable fact - ultimately rooted in the nonamenability of 
S0(1, N) - that this does not happen here. 



3.3 The structure of the ground state sector 

In fact the ground state sector of the S0(1, N) nonlinear sigma-models can be described 
very explicitly. We first summarize the result informally and then give a precise version 
in the form of a theorem. 

5*0(1, A^) nonlinear sigma-models defined on a finite d- dimensional spatial lattice have 
infinitely many generalized ground states transforming irreducihly according to the limit 
of the principal series. Every generalized ground state of the system lies in the linear 
hull of a single group orbit consisting of strictly positive functions. 

Note the sharp contrast to the ground state structure of the compact models: 

S0{N+1) nonlinear sigma-models defined on a finite d-dimensional spatial lattice have 
a unique ground state (which is a SO{N-\-l) singlet and which is strictly positive up to 



The precise form of the above statement is: 

Theorem 2: Let be the transfer operator (2.5), (2.6) of the S0(1,A^) nonlinear 
sigma-model. Then is a operator on — L^(]HI^^) with purely essential spectrum. 
There exists a unique function fio(^) with the following properties: it is strictly pos- 
itive and /^-invariant, i.e. p{A)Qo{n) = r2o(^)) for all A G S0(1,A^). Further Qo{n) 
is independent of one variable, Ux^ say, and square integrable in the other variables 
/ Hx^io'^^^^^) ^o(^)^ ^ terms of this function and V{n) := -ffo,o,o(^) ^ in 

(A. 10) the ground state sector ^(T^) of is given by 



In particular all generalized ground states Q G ^(T^) C L°° of transform according to 
the limit ttq of the principal series and are contained in the linear hull of a single group 



a phase). 




(3.30) 
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orbit. Explicitly, the former means they transform equivariantly according to 



n{A-^n) = (7ro(A)Q)(n) , A e S0(1, N) . (3.31) 

The theorem in particular guarantees the existence of generalized ground states in the 
sense defined in section 3.3. The crucial existence statement is that for the function 
Qo("^)- The proof of Theorem 2 is deferred to [37], where it appears as a special case 
of more general results. The application to the transfer operator of the nonlinear 
sigma-model rests on a technical Lemma which we present here: 

Let us denote by / : — a function satisfying fx{n) ■ fy{n) — rix ■ riy for all 
x,y E {1, . . . , Lg}"^ and fxoip) = for one Xq. Then 

/ H dn{nx) r^{n, fin); 2) < oo . (3.32) 

To see this we express the two-step transfer matrix in terms of the one-step Tg(n, n'; 1) 
defined in Eq. (2.4). This gives 

/"j]dO(n,)T^(n,/(n);2)= /"j] dO(n,)df7«) (3.33) 

X X 

X exp I - • {fix + fx{n)) + < • + ■ n^+f, - 4]| . 

X n^D 

We now estimate Ylx^xo'^^'x ' fx{n) > — 1 and view the result of the rix integrations 
as a function F[n'). It is p-invariant so that one of the spins can be frozen to a fixed 
value, say n'^^ — . The n'^^ integration then can be done and one obtains 

j J] dn{nx) Tp{n, f{n); 2) < L>^,^ e'^^^'-^) j J] dnin,) J] dn{n'^) x 

XyiXO 



X exp I - ■ < + < ■ n;+^ + ■ n^+f, - 3]} 



(3.34) 



Using '^^Tix-n'^ > rixQ-n^ + L'^ — 1 the integrals factorize into products of gauge-fixed 
one-dimensional partition functions and hence is manifestly finite. 

The property (3.32) guarantees that to our transfer operator (with only essential 
spectrum) one can associate a compact transfer operator (T^)o (with only discrete spec- 
trum) whose unique ground state wave function is the ^Q{n) featuring in Theorem 2; 
see [37] for details. 

As stated in the Theorem, the evolution operators of the nonlinear sigma-models given 
by (discrete Euclidean time) or Tr (continuous Euclidean time) both have purely 
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essential spectrum. The essential spectrum arises from the fact that or T,- cannot 
have generalized eigenstates transforming according to a finite dimensional irreducible 
representation of S0(1, N). The spectrum of these operators is a closed bounded sub- 
set of [0,00). Although there can be no normalizable eigenfunctions for the spectral 
values ||T^|| or ||Tt-||, it is not excluded that there exist infinite multiplets of normal- 
izable eigenfunctions corresponding to other spectral/eigenvalues. This situation would 
be interesting in that it might pave the way to a (quasi-) particle interpretation of the 
spectrum. In contrast to the universal structure of the ground state sector the existence 
or non-existence of such normalizable multiplets is a specific dynamical feature. In the 
case of the transfer matrix (3.16) the spectrum remains purely essential for any p in- 
variant potential; for certain potentials infinite multiplets of normalizable eigenfunctions 
might exist. In the absence of a potential, however, this is excluded. Indeed, from (3.14) 
we know that the kinetic part Tq and exp(— tHq), Hq := — ^Z^x^x'^' have absolutely 
continuous spectrum. 

The example of the kinetic part Tq of the transfer operator (to which all results of course 
apply in particular) can be used to get a feeling for 'how' the generalized ground states 
manage to be linear combinations of real and strictly positive functions: a complete 
system of real eigenfunctions is given by the tensor products of the functions (A. 10). 
Projecting out the SO^A^) singlet yields 

^oHy^,L{^) = Yi ^P^^i^^) Hui,L{n) , with 

H^,L{n) / d^{A) H H^^,i^,^^{An,) , (3.35) 

where the multi-index L refers to the set Ixjin^, x G A^, and uj_ = (cji, . . . ,ujlJ. On 
the other hand wc know from (3.7), (3.14) that the supremum of the spectral values 
of To is assumed if a;^; — > 0, for all x e Aj. The limiting eigenfunctions Hq l are 
real and strictly positive almost everywhere (a.e.). On the other hand an intcrtwiner 
Q{ui\uj) from Cn{(^i) ■ ■ ■ Cn{oJLs) to C^iuj) can be seen to contain 5{uj — Ylx^^:) 
as a factor. Assuming that Q{lo\uj) has a well-behaved limit, the irreducible component 
^a)(^) = (Q(^|'-^)-f^ui,L)(^)) will likewise be real and positive a.e. as — 0, in accordance 
with the above result. 

Theorem 2 is a special case of far more general results proven in [37]. Roughly speaking 
the above structure of the ground state sector turns out to arise mainly from the interplay 
between group theory and the general properties of a transfer operator. It thus admits a 
generalization largely independent of the details of the dynamics, which we outline here. 

Let T be any transfer operator in the sense defined in section 3.2. Then [37]: 

• T is a noncompact operator on with purely essential spectrum, (t(T) = (7ess(T). 

• Once the existence of a single strictly positive L°° ground state is guaranteed the 
ground state sector ^(T) of T assumes a certain universal form, independent of 
the details of the dynamics! 
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There exists a transfer operator Tq, uniquely associated with T and with the 
same spectral radius, such that the ground state sector of G (T) is related to that 
of g{To) by 

^?(T) ~Span{QoWP(^n,J, AeS0(l,7V), e ^?(To)} • (3.36) 



• The (generalized) ground states of Tq are singlets, so that by (3.36) all generalized 
ground states of T transform according to the limit ttq of the principal series: 
n{A-^n) = (7ro(^)n)(n) , A e S0(1, A^). 



3.4 Thermodynamic limit, SSB, and time-slice be 

In a hamiltonian formulation the thermodynamic limit is hard to control because the 
'Hilbert space changes'. The way to proceed is to take the thermodynamic limit on the 
level of the correlation functions and then reconstruct a Hilbert space formulation via 
a Osterwalder-Schrader reconstruction. We return to the second aspect in section 3.5. 
Of course even on the level of expectation values the thermodynamic limit is difficult to 
control. Interestingly there is an elegant argument saying that the limit of the functional 
measures (whenever it exists as a mean) cannot be invariant for all D > 1. 

Theorem 3: Expectations ( )A,/3,bc of the SO(l,iV) nonlinear sigma-model defined on 
a finite D-dimensional lattice A with SO^(iV) invariant be and gauge fixing cannot have 
an S0{1,N) invariant thermodynamic limit ( )oo,/3,bc '■— limA_»zo( )A,/3,bc- Specifically, 
there exist bounded continuous functions 0{n) of one spin such that 

(C»(^n))oo,/3,bc 7^ (0(n))oo,/3,bc for some ^eS0(l,7V). (3.37) 



As noted in section 2.3 for D > 3 a similar conclusion was reached in [21] by very 
different means and based on a different criterion. In the present setting the symmetry 
breaking is essentially a consequence of the fact that a nonamenable group does not 
have an invariant mean over 'nice' function spaces [40]. Here we consider the space of 
bounded continuous functions C6(S0(1, N)) on the group manifold. Equipped with the 
sup- norm it forms a commutative C*-algebra with unit, so that the usual concept of a 
state applies. 

The proof of Theorem 3 is a straightforward generalization of the argument origi- 
nally presented for = 1 in [18]. The expectation of a single-spin observable at 
X = {xi, . . . jXd) can be written as 

{0{n^))A,i3,hc = / diJ,A^bc{n;x)0{n^) . (3.38) 
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By construction, for all finite lattices A the one-spin measure (i//.A,bc(^; a;) , x e Aa;^, is 
a normalized probability measure depending parametrically on L^. One can thus view 
(3.38) as bounded, positive, and normalized linear functionals ('states') on the func- 
tions in Cb{^0{l,N)) which happen to be independent of the variables in the SO^(A^) 
subgroup. By the theorem of Banach-Alaoglu [38] there is therefore a subsequence 
of lattices A on which the states ( )A,/3,bc converge to a limiting state ( )oo,/3,bc on 
Cfc(SO(l, iV)). Because S0(1,A'") is non-amenable this limiting state cannot be invari- 
ant. There must exist functions Q G Cb{^0{l,N)) such that their average in the state 
( )oo,/3,bc is noninvariant. For all finite lattices {Q{B))A,i3,hc = (C^(^^))A,/3.bc holds, where 
0{ni n eMn = S0{1,N)/S0\N), is the SO^(iV) average of the function Q{B) and 
n = B SO^(A^). Since both the observable and the sequence of states are SO^(A^) invari- 
ant, the limit will also be invariant, (Q(-B))oo,/3,bc = (C'(n))oo,/3,bc- On the other hand by 
definition of Q one has (0(An))oo,/3,bc = (<5(^-B))oo,Abc 7^ (<5(-B))oo,/3,bc = (C(?T'))oo,/3,bc, 
for some A e S0(1, N), as claimed. 

We remark that when the continuity requirement on the symmetry breaking observable 
is dropped, the key step in the argument follows more directly from the known character- 
izations of an amenable symmetric space. A symmetric space G/ H (G a locally compact 
group and H a maximal subgroup) is called amenable if there exists a G-invariant mean 
on L°°{G/H) (see [41]). A unitary representation tt of a locally compact group G on a 
Hilbert space 7i is called amenable in the sense of Bekka if there exists a positive linear 
functional u over B{T-C) (the C*-algebra of bounded linear operators on Ti.) such that 
u{'ri{g)T'ri{g)^^) = uj{T) for all g E G and all T G B{l-i). Then the following three state- 
ments are equivalent: See [42] and [43] (i) G/H is an amenable symmetric space, (ii) 
the quasiregular representation pi of G on L'^{G/H) is amenable in the sense of Bekka. 
(ill) The quasiregular representation pi almost has invariant vectors in the sense that 
for all compact K G G and all e > there exists a unit vector ip G L'^{G/H) such that 
||pi((?)'?/^ — t/'II < e. Applied to the case at hand we know that the quasiregular representa- 
tion pi of S0(1, A^) on Hat does not almost have invariant vectors, e.g. from the explicit 
decomposition (3.4). Thus H^r is a nonamenable symmetric space and by the above 
argument there must be symmetry breaking observables O G L°°(E[jv), i.e. essentially 
bounded and measurable functions O of one spin such that (3.37) holds. 

In the rest of this section we present a nonrigorous argument that these symmetry 
breaking single spin observables are the rule rather than the exception. To this end we 
introduce expectations with a third type of boundary conditions which take advantage 
of the extremal configurations G H^" . As stated after Theorem 1 the existence of 
these extremal configurations, although unproven at present, is highly plausible for the 
transfer operator of the noncompact sigma-models. The definition of the expectations 
( )a,/3,3 based on these configurations is as follows. We set 




(3.39) 
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The expectations referring to n* are defined for a one-spin observable by 



J Ip[n^,n^,Lt) 

for a two-spin observable by 

(0K,n,))-^,3 = — - / n n Mny)x (3.41) 

ip[n„n„i.t) J ^^^^^ ^^^^^ 

Tp{n^, n^; Lt/2 + xd) 0{n^, Uy) Tp{n^, Uy] yo - XD)Tp{ny, n*; Lt/2 - yn) , 

and so on. The notation in (3.40), (3.41) is the same as in (2.13); compared to the second 
type of be the lattice now ranges over time slices A^;^ ,xd = —Lt/2, . . . ,0, . . . ,Lt/2, 
with Lt even; all spins in the time slices A±l^/2 are frozen to the special configuration 
n*. Although we expect to be SO^(A^) invariant (in that all spins can be chosen to be 

equal to n^) the be also work arbitrary e The invariance under S0\N) (which 
was a feature of the other two types of be) then has to be restored by group averaging. 

The advantage of these boundary conditions is that by Theorems 1 and 2 the Lf ^ oo 
limit can be analyzed in a similar way as in the 1-dimensional model [18]. For example 
for S0(1, N) invariant observables one has 

lim (0(n,,n,))^^3=||T|r^-^^[27r^/2r^^/2)f]o(n,)]-i / TT dn{n,) x 

Lt— >oo I 



xO{n^,n^)Tp{n^,n^;yD - XD)^o{na:)V{{n^)xo ■ n^o) ■ (3.42) 

Here the limit is taken on a subsequence {Lt)j^^ C N as in Theorem 1, and Jlo(^) 
is the p- invariant positive function in Theorem 2. The normalization is such that the 
limit functional obeys limLj^oo(l)A*/3 3 ~ ^- This line of argument clearly generalizes 
to the expectations of observables depending on any finite number of spins. For such 
observables also the subsequent Lg ^ oo limit exists on subsequences, by the theorem of 
Banach-Alaoglu. Clearly this argument does not depend on the number of dimensions. 
We conclude: 

The expectations (C)a*/3 3 of all local 5*0(1, N) invariant observables, defined with n = 
boundary conditions at x^ = '^Lt/2, have a pointwise finite and explicitly computable 
thermodynamic limit, \imL^^oo^^^Lt^oo{0)^*p ^. 

For noninvariant observables the evaluation of the thermodynamic limit is more difficult. 
An exception are observables depending on a single spin only. We shall now argue 
that basically every nontrivial bounded function of a single spin will signal spontaneous 
symmetry breaking in the sense of (3.37). This can be seen when using type 3 be 
combined with a slightly heuristic use of Theorems 1 and 2. To this end consider the 
family of measures d^^^^{n]x) in (3.40) with type 3 be. Let us write Tp^n^^n^] L^) — 
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0{d{Lt)) for the leading asymptotics of the denominator. The density of the measures 
diiA,3{n;x), X G A^.^, then behaves as 

d{Lt)no{nfV{{n,)^,-n^,f for ^ oo . (3.43) 

On general grounds d{Lt) ^ as oo [37]. The density (3.43) thus vanishes 

pointwise in the limit. The proof of Theorem 2 is based on the fact that the state space 
and the unitary representation p can be factorized according to [37] 

~ L'^{mN)®L'^{^f), with M = / p{SO{l, N)) , 
p ~ Pi^pL, (3.44) 

Here L^(H[jv) is the factor carrying the dependence on the preferred variable n^^Q, and 
L'^{J\f) is the factor invariant under the action of p. We identify p with pf^% the Lg- 
fold tensor product of the quasiregular representation and p^^^ is a representation acting 
trivially on Li^^M). In particular VLq lies in L'^{J\f) and thus has finite norm ||f2o||A/' with 
respect to the invariant part of the measure. Based on this factorization and (3.40), 
(3.43) one obtains 

lim ((!:'(n^))A,^,3 = ||Qo||^ lim j dpLt{nxo) O{n^o) . (3.45) 

Here dpL^iuxo) is a one-spin measure whose density scales like d{Lt) V{{n^)xo ■ n^oY for 
Lt oo. The point here is that the second factor on the right hand side of (3.45) is 
independent of Ls and, whenever the limit — oo exists, it is not S0(1, N) invariant. 
Assuming that this is the case one can take the Lg ^ oo limit of Eq. (3.45). This 
affects only the Hl^ollAr term which is p invariant for any finite Lg and hence also in the 
limit. The second factor however is noninvariant and gives rise to spontaneous symmetry 
breaking - practically for every nontrival bounded one-spin observable, as asserted. 



3.5 Osterwalder-Schrader reconstruction 

The purpose of the Osterwalder-Schrader reconstruction is to reconstruct a Hilbert space 
and transfer operator as well as a translation invariant state ('vacuum') Q from correla- 
tion functions (or more generally expectation values) in the thermodynamic limit. Since 
the infinite volume limit for the transfer matrix cannot be taken ('the Hilbert space 
changes'), this is the only way in which a physical interpretation of the model in infinite 
volume can be achieved. The general procedure has been described in many places, see 
for instance [44, 45, 46]. The following discussion applies to S0(1,A^) sigma-models in 
any dimension D > 1. For the 1-dimensional version the Osterwalder-Schrader recon- 
struction is discussed in detail in [18]. 

The crucial properties required in a lattice system are 
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(RP) Reflection positivity, and 
(TI) Time translation invariance. 

The property of RP in our model can be stated as follows: Denote by C+ a 'suitable' 
linear space of continuous functions of finitely many spins rix at positive 'times'. If 
O G C+, let i^C be the complex conjugate of the same function of the time reflected 
spins. Then 

(WO)A,^,bc > . (3.46) 

It is satisfled by our model as long as the volume is flnite, the temporal size Lt is 
even and we use boundary conditions that are time-symmetric; this follows from the 
representation of the system in terms of the transfer matrix T, see Section 2.1. For 
example it holds in the fixed spin gauge (be = 2) with periodic be in time direction 
if the fixed spin is chosen to have time coordinate Lt/2 (identified with —Lt/2). It 
also holds for the fixed time-slice gauge (be = 3) considered in (3.41), (3.42). In the 
latter case the existence of a thermodynamic hmit is guaranteed at least for S0(1, A?") 
invariant observables. Translation invariance is manifest already on a finite lattice and 
thus holds trivially also in the limit. Since for S0(1, N) invariant functions the different 
gauge fixes are presumed equivalent this should entail the existence of a thermodynamic 
limit also for i = 1, 2 be, where the limit is then given by the same formulas as for the 
type 3 be. For the fixed spin gauge (i=2) spatially periodic be are used; translation 
invariance is not manifest on a finite lattice but should be restored in the limit. We 
remark that for SO(l,iV) noninvariant observables, such as our 'Tanh' order parameter 
the properties (RP) and (TI) are not obvious, since in on a finite lattice (RP) is violated 
in the translation invariant gauge, while (TI) is violated in the fixed spin gauge and in the 
fixed time slice gauge. Nevertheless it is reasonable to assume that both properties are 
restored in the thermodynamic limit. For noninvariant observables in C-f we therefore 
assume here the existence of a translation invariant thermodynamic limit. We write 
( ■ ) for the hmiting functional ( ■ )oo,i3,h ^ = 1)2,3, with the above specifications. By 
construction it then also has the property RP. 

The property (RP) allows one to construct a Hilbert space both on a finite lattice and 
in the thermodynamic limit, which we denote by TYa and Tios, respectively. The defini- 
tions {0',0)a := (C^?0')A,/3,bc and {0',0)os ■= {O^O') define a positive semidefinite 
scalar product on C+. By dividing out the subspace of elements with vanishing norm 
{0,0)x,p,hc = and {0^0)os = Oi respectively, and completion we obtain the Hilbert 
spaces and Ti-os, as described. Importantly Ti-os will in general not be separable. 
This was found explicitly in the 1-dimensional model, and it is unlikely that separability 
will be restored by 'adding' spatial dimensions. For example the OS reconstruction of 
the solvable noncompact model of a massless free field in two dimensions likewise leads to 
a nonscparablc state space [47]. In contrast the spaces Ti.\ are, for any finite |A| = LgLt, 
isometric to L^. Denoting the isometry by V\ : TYa ^ L^, the unitary representation 
p on induces one on TCa, namely pa := ^V^a- Our second assumption is that 
the thermodynamic limit of Pa{A) exists weakly, i.e. the limit lim|A|^oo(PA(^)C', C'Oa 
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exists for all O, C G C+, and defines a measurable function of A G SO(l,iV). This 
defines a measurable action pos of S0(1, N) on Tios- Guided by the properties of the 
1-dimensional case, we do not expect or require this action to be continuous. Further 
Pos is expected to be unitary only on a closed subspace TVos of T^os- Since Tios is 
in general not separable an alternative described by Segal and Kunze applies; see [48] 
and [18] in the present context. The upshot is that 'Hos decomposes into a direct sum 
1~Cqs — '^os ® ^05' where the restriction of pos to Ti-os '^os i^ continuous and 
singular, respectively. Here singular means that {ips, Pos{^)'4^s)os — for almost all 
A G S0(1, A^) and all ips G '^os- "^os i^ separable, Ti^s i^ absent. In one dimension 
such a representation pos could be constructed explicitly, and Tipg turned out to be 
nontrivial. The explicit form of pos also entailed that the fl induced by 1 G C+ is 
actually an element of a ground state orbit {pos{A.)^l, A G SO(l,iV)}. The infinite 
dimensional closed subspace of Tios spanned by this orbit was contained in Hq^, that 
is, the action was continuous and unitary. 

The nonamenability of S0{1,N) now has no direct bearing on the existence of almost 
invariant vectors for pos, since even a nonamenable group can have amenable represen- 
tations. However unitary amenable representations vr (in the sense of Bekka, defined in 
the remark following Theorem 3) are characterized by the fact that 7r®7f almost has in- 
variant vectors, where tt is the conjugate representation ([42], Theorem 5.1). We expect 
that this can be used to rule out that pos restricted to Ti-os i^ amenable. Schematically, 
we identify pos with a weak limit pos = w — lim|A|^ooPA (provided it exists) of the 
unitary representations p\ on Ha as described above. Then, using — p — pi® p^v 
from (3.44), one formally has 

Pos ® Pos ^Pi®pi®(w- hm p^^ (8) p^\ . (3.47) 

\ |A|^oo / 

Since by assumption pos defines a measurable representation of S0(1, N) on Hos and 

Pi®Pi does so on L^(EI^), the second factor in (3.47) should likewise define a measurable 
representation of S0(1, A^) on some (nonseparable) Hilbert space. Since pi ® pi does not 
have almost invariant vectors it is highly plausible that pos restricted to V-os cannot 
have almost invariant vectors. 

Next we address the reconstruction of a transfer operator. Translation r by one lattice 
unit in positive time direction maps C+ into itself. Since by assumption the limiting 
expections ( ) are translation invariant it can be shown by standard arguments (see 
[44, 45, 46, 18]) that this map lifts to a well-defined bounded symmetric operator Tos 
on Tios-i which is the desired reconstructed transfer operator. By construction, there 
is a normalizable state Q induced by the constant function 1 G C_|_, which is a proper 
eigenstate with eigenvalue unity of Tos ~ in sharp contrast to T which had no proper 
eigenstates in ~ Ha. We denote by QiTos) the set and closed hnear subspace of all 
(normalizable) ground states of Tos in Ti-os- 

Concering the interplay of Tos with pos there are two main cases to consider: First the 
weak limit defining pos does not exist even when restricted to Q{Tos)- Second, pos does 
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exist at least when restricted to Q{Tos), in which case the action of pos can be unitary or 
nonunitary. In the 1-dimensional model the second possibility was realized with a unitary 
action. Indeed Tqs commuted with pos on all of Tlos- Further Q{Tos) was a subspace of 
'H.Qg and pos restricted to Q{Tos) was equivalent to ttq, the limit of the principal series. 
The first possibility should be taken into account based on experience with spontaneous 
breaking of compact symmetries in dimensions D > 3. In this case the group no longer 
acts on the (unique) ground state because of the infinitely many degrees of freedom that 
would have to be transformed. In the case of a compact symmetry one has the option 
of averaging the expectation values over the symmetry group, thereby introducing a 
'large' (but still separable) Hilbert space as a direct integral over the pure phases. In 
this Hilbert space there is then degeneracy of the vacuum and the symmetry group 
acts nontrivially on the vacuum space. The original vacua are recovered by an ergodic 
decomposition of the symmetric state. Because of the non-amenabihty of S0(1, N) one 
does not have this option here. In summary, we can envisage the following scenarios for 
the interplay between pos and Tqs' 

(1) Pos does not exist even when restricted to Q{Tos), or on the restriction Tqs and 
Pos do not commute. 

(2) Pos does exist at least when restricted to Q{Tos), and on this subspace Tqs 
and Pos commute. Q{Tos) then decomposes into an orthogonal sum of subspaces, 
Q'"'(Tos) on which pos acts unitarily and Q"'^(Tos) on which pos acts nonunitarily. 
Further pos restricted to Q'"-{Tos) is expected to be nonamenable and Q^{Tos) 
decomposes in to a direct sum G'^{Tos) © G'^i'^os), where the restriction of pos is 
continuous on Q'^{Tos) and singular on S^{Tos)- One or both of these subspaces 
could be trivial. 

(2a) If Q'^{Tos) is nontrivial it carries a unitary continuous representation of S0(1, N), 
which one can assume to be irreducible. Based on the results of Section 3 a 
plausible candidate is again the limit of the principal series ttq. If Q'^i^os) is 
nontrivial the group acts on it discontinuously as a 'permutation group'. Such an 
exotic situation was found in the ID case for a certain non-vacuum subspace of 
'Hos- 

At present we do not have enough information to determine which of the above scenarios 
holds. All however represent refinements of the fact that the symmetry is spontaneously 
broken. 
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4. D = 2 : Numerical simulations 



Although well suited to address structural issues, the hamiltonian formalism used in 
sections 2 and 3 is not ideal to obtain quantitative results. While in the compact models 
this is still feasible [49, 50] the intricate group theory required in the noncompact case 
seems to render such an approach unattractive for models with a noncompact symmetry. 
In Sections 4 and 5 we therefore study the dynamics in terms of correlations functions, 
first by numerical simulation and then via the large N expansion. 

We have performed simulations of the S0(l,2) sigma-model on square lattices of linear 
dimension L = Lg = ranging from L=20 to L=128. The simulations were performed 
at different coupling values, and with the two choices for the gauge-fixing described 
in Section 2.2. We now describe briefly the Monte Carlo algorithms employed in the 
simulation of these averages (C)A,/3,i, i = 1, 2: 

1. For the average (0)a,/3,i (translationally invariant gauge and periodic be) a Monte 
Carlo sweep through a L x L lattice is defined as Metropolis updates of randomly 
chosen pairs of spins , according to 

fixi ^ fixi + r, fix^ ^ fix^ - r , (4.1) 

where the two-dimensional vector r = r(cos 0, sin 0), r.cf) chosen randomly in the 
ranges (0, rinax)i (0, 27r) respectively, and with r^ax adjusted to yield acceptance 
rates close to 50%. The symmetric update of pairs of spins ensures that the gauge 
constraint Ylix'^^ = is preserved. The initial configuration (we typically take 
cold starts, with fix — 0) is of course chosen also to satisfy this constraint. The 
proposed update (4.1) is then accepted or rejected on the basis of the change in 
the effective action 

,Si = ^ ■ Ux+f, - 2 In 5^ n° + ^ In n° . (4.2) 

After an initial run of 100000 sweeps to equilibrate the system, configurations 
are stored subsequently at intervals of 2000 sweeps, which is than adequate to 
decorrelate the configuration (for the observables measured, typical autocorrelation 
times are at most a few hundred sweeps). The results presented in this paper derive 
from averages over ensembles of 5000 independent configurations, unless otherwise 
stated. 

2. For the averages (0)a,/3,2) a sweep is defined as updates of a randomly chosen 
single spin n^.^, not including the fixed spin at Xq, 

^xi rixi + r, xiy^xo, (4.3) 
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with the random vector r chosen as above, and the effective action in this case 



Again, we typically followed an initial equilibration with 100000 sweeps by 5000 
measurements spaced at 2000 sweep intervals. 

We have found that long simulations of this model lead to unrehable results unless a 
high quality random number generator is employed. Specifically, violations of translation 
invariance of two-point functions at the 3 to 4 standard deviation level were found when 
the random() function packaged with GNU gcc (specifically, gcc-2.96) was used. The 
ranlux generator developed by M. Liischer [51], double precision, and with the luxury 
level set to 2, was used in all simulations reported in this paper. With this generator, 
we have found no violations outside of statistics in expected symmetries of measured 
observables. 




Figure 1: Typical configuration of the 'heights' n^, x G A, at strong coupling (3 = 0.1 for 
L = 64. Blue, green, orange corresponding to low, medium, high values of n^, respectively. 
The mean value is (nP) = 5.11. 

It is instructive to look at some typical configurations in the parametrization Ux = 
(^x, X G A. As discussed in Section 2.2 with the translation invariant gauge 

fixing one expects the S0^(2) subgroup to be unbroken. The compact spins will 
then be distributed similar as in the massless phase of the familiar 0(2) model. The 




(4.4) 
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Figure 2: Typical configuration of tlie 'lieights' n^, x G A, at weak coupling /3 = 10 for L = 64. 
Blue, green, orange corresponding to low, medium, high values of n^, respectively. The mean 
value is = 1.067. 

novel feature are the noncompact components C,x = for which we show some typical 
configurations at weak and strong coupling in Figs. 1,2. One sees that at strong coupling 
the mean value (^i°)a,/3,i is large, with relatively large localized fluctuations rendering 
nearby spins almost uncorrelated. For weak coupling on the other hand most of the 
spins are 'frozen' close to bottom of the hyperboloid, (n°)A,/3,i ~ 1, and nearby spins are 
correlated, both in 'height' and in direction n. 



4.1 Spin two-point function and energy correlator 

The spin two-point function {n^ ■ ny)\^p^i is the simplest S0(1, 2) invariant bilocal object 
constructible in the model. The thermodynamic limit of this quantity can be studied 
numerically by simulating various size lattices at flxed f3. The results at /3 =10 for square 
lattices of linear size L = \/\A\ = 32, 64, 128 and i = 1 (periodic boundary conditions, 
translationally invariant gauge) are shown in Fig. 3. They suggest the existence of a 
flnite thermodynamic limit, consistent with the analytical arguments in Section 5. It also 
illustrates that the spin two-point function increases with increasing separation \x — y\. 
This somewhat peculiar behavior has been observed before [2] and can be understood 
analytically both in the ID model [18] and in a large N analysis, see Section 5. 

Another natural invariant observable is the 'energy' or action density = 2(1 — ■ 
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Figure 3: Two-point function (n^ • ny)\^p^i, for /? = 10 and varying L. 



i^x+ji)- Its expectation value appears in the invariant combination of the Ward identities 
(2.22). Here we study the connected part of its two-point function and probe for nontriv- 
iahty and clustering. The subtractions involved in extracting the connected part involve 
large cancellations, and we have had to perform very long runs (collecting ensembles of 
40000 configurations) on somewhat smaller (20 x 20) lattices to find a meaningful signal. 
The connected part was also very small in the weak coupling regime, so we needed to 
go to strong coupling; the results of Fig. 4 correspond to /3=0.1. At least for separa- 
tions r < lattice spacings there is then a nonvanishing signal. The fact that the 
signal disappears so rapidly makes it impossible to draw any firm conclusions from the 
numerical data on the nature of the asymptotic falloff: for example, to distinguish be- 
tween the power behavior suggested by naive dimensional reasoning, or exponential 
falloff. In summary, we find a nontrivial energy correlator rapidly decreasing at nonzero 
separations. 



4.2 Two-point function of the Noether current 

The Ward identity (2.27) for the longitudinal momentum space current correlators pro- 
vides a stringent test that the simulation scheme is fully respecting the symmetries of 
the model. In Fig. 5 we show the comparison of the left and right hand sides of (2.27) 
on lattices of size 32 x 32 and 64 x 64 (periodic boundary conditions, translationally 
invariant gauge), for the boost {{ah) =(01)) and rotation {{ah) =(12)) Noether currents. 
The agreement is within statistical errors, except for the lowest momentum modes. In 
fact, we show in Appendix B that at fixed nonzero momentum the delta function gauge 
constraint induces a finite volume correction of order 0{h\V/V) to (2.27). These finite 
volume corrections are largest at the edge of the Brillouin zone, i.e for the momentum 
modes of order p ~ 1/L (see Appendix B, Fig. 12). The transverse Noether correlators 



37 




1 2 3 4 5 

r=|x-y| 

Figure 4: Connected correlator {E^Ey) - {E^){Ey), 20 x 20 lattice, (3=0.1. 

are nontrivial, and are shown in coordinate space in Fig. 6. The falloff is roughly as 
expected on dimensional grounds (fits also indicate a logarithmic component ln(//r)/r^). 



4.3 Tanh observable — spontaneous symmetry breaking 

Finally, we have used our generated ensembles to measure the order parameter T(g) 
introduced in Section 2.4 as a signal of spontaneous symmetry breaking of the SO (2, 1) 
group. For N = 2 the average Tq{^) in Eq. (2.25) is a strictly decreasing positive function 
for ^ e [1, oo), with the limiting values Tq{l) = tanh{^/q^^^^) and Tq{oo) as below. As 
described in section 2.3 by a convexity argument one expects 

(Te(n))A,^,i > T,((n°) A,;3,i) > Tg(oo) = 1 - - arccos ^1 - q'^ , (4.5) 

TT 

where the limit is taken from [18]. In Figs. 7, 8, the results for the three functions in 
(4.5) are shown for weak and strong couplings, respectively. One sees that for weak 
coupling the spins are almost 'frozen' and {Tf,{n))A,/3,i practically coincides with the 
average Tq((n°)A,/3,i)- For strong coupling, on the other hand, genuine dynamics sets 
in and both quantities differ. Importantly, in cither situation the symmetry breaking 
is manifest in that T(g) := (Te(n))A,/3,i is a nontrivial function of q. Tq{l) vanishes on 
account of the unbroken S0^(2) symmetry; the lower bound Tq{oo) guarantees that the 
curve cannot 'flatten out' and vanish identically in the thermodynamic limit. 

As discussed in Section 2.4 the divergence of (n°)A,/3,i as |A| — > oo is not really a test of 
spontaneous symmetry breaking. Moreover, because of the expected soft (logarithmic) 
divergence a very large range of lattice sizes would be needed in order to pin-down 
a suspected divergence of (n'')A,/3,i. For example at /5 = 10 we obtained (n°)A,/3,i = 
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Figure 5: Ward identity for longitudinal Noether currents: L = 64, /3 = 0.1 (top) and /? = 10 
(bottom) . 



1.0585,1.0695,1.080, for L = v^|A| = 32,64,128, respectively. In itself this would 
hardly constitute convincing evidence for a divergence in the L — > oo limit. 
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Figure 6: Transverse Noether current correlator, (3 =10, L =32,64. 



5. D = 2 : Large N analysis 



The noncompact SO(l,iV) sigma-models may be analyzed in the usual large limit, 
i.e. N ^ oo, X := N/(3 = g^N fixed, by saddle-point techniques analogous to those 
used in the compact case. However, several important differences arise which alter 
quahtatively the results in the noncompact case. The large N analysis is especially 
useful for examining qualitative features like the behavior of correlation functions in 
the thermodynamic limit, thereby providing a guideline for the correct extrapolation 
of numerical results to the continuum limit. We adopt the setting of Section 2.1 and 
perform a large N analysis of the lattice-regularized model with periodic boundary 
conditions, using the translationally invariant gauge- fixing 1 in Section 5.1 and the 
fixed-spin gauge 2 in Section 5.2. We consider square lattices only and write ( • )L,N/x,i 
for ( • )A,/3,i, ^ = 1,2, with L — ^/\A'\ the hnear size of the lattice. 
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Figure 7: Order parameter T{q) for weak coupling. 
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Figure 8: Order parameter T(q) for strong coupling. 

5.1 Lcirge N analysis in a translationally invariant gauge 

By (2.10) the partition function Zi = Zi{A,p = N/X) has the form 

Zi ^ Jl[dn,5inl-l)5(^J2'^-) 

X X 

Implementing the nonlinear constraint as usual with an auxiliary field a^i (5.1) becomes 

X ^ ^ x,n XX ^ 

X j Y{dn^^^v\- n(^Y^ n^{-A)^yny + i ^ a^n^ j I , (5.2) 
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with Aa;y the discrete lattice Laplacian. On integrating out the n-field one finds, up to 
irrelevant multiplicative factors 

Zi ~ / dn^da^ exp {-NSi{n^, a)} , (5.3) 
with the effective large N action 
Si^-^J2 - 1^ E < - ^ E - 1] + ^Tr' In [-A + 2tXa] . (5.4) 

xy XX 

The prime on the trace in (5.4) denotes omission of the zero mode of the Laplacian, in 
keeping with the Si^^fix) constraint in (4.1). 

Note that in contradistinction to the compact case, the n° field is not integrated out in 
defining a large effective action, and we are led to the problem of determining a joint 
saddle-point in (n°, a) field space. We now show that there always exists at least one 
translationally invariant joint saddle-point of (5.3). Let 

nl = n + ifjj:, ax = -ia + , (5.5) 

where n, a are real, but the phase of the fluctuation variables rjx, Cx (despite the sugges- 
tive notation) remains to be determined later by a detailed analysis of the local structure 
of the saddle-point. The auxiliary field integrations in (5.3) run initially along the real ax 
axes but (as is frequently the case in Hubbard-Stratonovich type saddle-points [52, 53]) 
require a deformation through a purely imaginary saddle-point. In contrast, the rem- 
nant field integrations over n° run along the semiaxis [l,oo), and we shall find real 
saddle-point (s) with n > 1 and real. 

The saddle-point conditions — §^ — ^ yield 

n'^l + X{-A + 2Xa);^, a = (5.6) 

where V := |A| = is the lattice volume. Note that a < 0, corresponding to a negative 
dynamically generated squared mass in the gap equation (5.6). Explicitly this becomes 

^2^^(l-cosp^)-^ Vn^ 

The discrete lattice momenta arc = ^^n^, = 0, 1,...,L — 1. Due to the infrared 
divergence of the sum in (5.7) the expectation value of the nP field diverges logarith- 
mically for V — > oo, specifically as {{n^)^) ~ ^logl^. Thus the dynamically 
generated negative squared mass in the gap equation is actually of order in the 
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N 


(nO) 


20 


3.980 ± 0.007 


40 


3.808 ± 0.004 


80 


3.651 d= 0.002 


160 


3.533 ± 0.002 


320 


3.458 ± 0.001 


640 


3.437 ± 0.0005 


OC' (sa(l(llo-i)()iiit) 


3.1136 



Table 1: Comparison of {n^)L,N/x,i for L = \ = 2{) from simulations with large N result. 

thermodynamic limit. Solutions of (5.7) with f[z) > correspond to n > 1. For any 
y, A there is always a root with 2; < 4 sin (^)^ and f{z) > 0. In the weak couphng regime 
defined by the inequality A < 4L^ sin (f 40 for large L) it is easy to see that this is 
the only root yielding n > 1. Henceforth we shall assume weak coupling, in the sense of 
the above stated inequality, and dominance of the single saddle-point with n > 1. 

We have performed explicit numerical simulations of the S0(1, N) model at values of 
ranging from = 20 up to = 640 on a 20 x 20 lattice to check the saddle-point result 
(5.7). The coupling was chosen in the weak coupling regime in the sense indicated above, 
specifically A =20. The convergence to the large N hmit is shown in Table 1. A fit of 
the data to a functional dependence of the form A + B/N + C/N'^ gives A =3.402. The 
numerical evidence suggests that the large N functional integral is indeed dominated by 
the saddle-point located in (5.6). 

One may also study field correlators in the large limit (assuming again the dominance 
of the saddle-point exhibited above). For example, the S0(1,A^) invariant two-point 
function to leading order in large A'^ takes the simple form 

{nx ■ ny)L,N/\,i "^n^ - XD^y , (5.8) 

with D := [—A -|- 2\a]~^. Of course {ux ■ nx)L,N/x,i = 1 as a consequence of the gap 
equation (5.7). The first term in (5.8) arises from the correlator, while the second 
term represents the cumulative effect of A^ 'spatial' n correlators, each of order 1/A^, with 
a negative sign from the indefinite metric. As the ft correlators fall off with distance, 
the S0(1, N) correlator evidently is an increasing function of separation. In Fig. 1 we 
compare the saddle-point result (5.8) for this correlator with simulation results obtained 
at finite A^ on a 20 x 20 lattice. 

Although we shall not compute 1/A^ corrections to the leading large A^ results here, it 
is of interest to study the local structure of the saddle-point in this theory, which is a 
prerequisite to evaluating the Gaussian fiuctuation corrections. Indeed, the routing of 
the field integrations through a joint saddle-point of the type studied here is somewhat 
more intricate than usual: we shall find that relative to the implied contour deformations 
in (5.5), further contour rotations are required in nonlocally defined field components. 
Expanding the large N action S{n + irjx, —ia + ^x) to second order in the fiuctuation 
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Figure 9: Large N convergence of invariant 2-point function, L = A = 20. 



fields "qxiixi one finds the quadratic form 

xy X X xy 

(5.9) 

with M := —A + 2\a. The appropriate routing of the field integrations through the 
saddle (5.5) is best analyzed by going over to momentum space: we replace J Yl^ dVxd^x 
by J Y[pdvip)d^{p) where rj^ := yJ2pV{p)^^^ '^^ etc. The quadratic form (5.9) now 
becomes 

5'='42:(I^^WI,WI^ + 2^,W5(-p))-^,(0f + ^ E 

(5.10) 

with M(p) :— 4 sin (^) +2Xa. We shall henceforth neglect the term involving 77(0)^ 
in (5.10), as it is of order 1/V relative to the rest. Defining the one-loop polarization 
function 

the quadratic action can be written as 

In (5.12) we change field variables from (?](p),C(p)) to (^(p) •= vip) + i^^C(p)) C(p))- 
This functional change of variable has unit Jacobian but is of course highly nonlocal in 

coordinate space. The quadratic action can now be expressed in terms of the variables 
57r(p) := (57 (p) +fi{—p))/2,fij{p) := {fi{p) — fi{—p))/2i, and similarly for C,- The integrals 
over these variables run initially along the real axis, but must be rotated in passing 
through the saddle-point depending on the sign of M{p) and ]\|^^^ — A^n(p) in (5.12). 



44 



For p =0 one has M(0) < 0, n(0) > 0. Accordingly the zero mode of the fj field must be 
rotated by 7r/2 in passing through the saddle-point, while the zero mode of the ^ field 
retains the initial contour orientation. The situation is reversed for nonzero momenta: 
namely, we find 

2A 

M{p) > , j^n^ > X^U{p) , for p 7^ , (5.13) 

which implies that for nonzero momentum modes the ^ contours should be rotated by 
7r/2, while the rj routing is unchanged. Calculations to next to leading order in 1/N 
necessarily require careful attention to the phases induced by these contour rotations. 
We again emphasize that the above statements hold for the case of a single dominant 
saddle point, when A < 4L^sin(^)^. 



5.2 Large N analysis in a fixed-spin gauge 
By (2.11) the partition function Z2 = Z2(A,/9 = N/X) reads 

^2 = / II (^^x^i^l - 1) ^i^xo - ^^) 

X exp J](n°+, - - ^ J](n.+^ - n^)A \ . (5.14) 

^ ^ X,IX X,fJ, ' ^ 

Introducing the auxiliary field a as usual, the partition function becomes in this case 

Z2 = y JJ (in°(ia^exp{-Ar52(n°,a)}, (5.15) 

with the effective large N action 

-52 = E n°(-A).,n; - r J] a.[(n°)2 - 1] + ixr'ln [-A + 2.Aa] , (5.16) 

xy X 

where in (5.16) the prime on the trace now implies a projection corresponding to omission 
of the integral over the field variable n^^. The presence of a fixed spin at a definite point 
on the lattice now means that the saddle-point will involve (n°, a) fields with a nontrivial 
spatial dependence. Writing in analogy to (5.5) 



nl = n^ + iVx: o:^ = -ioix + L: (5-17) 
ds _ as 



the saddle-point conditions ^ — 4^ — now lead to 



nl = l + XD^^, = -i-^(An)^, (5.18) 

zA Ux 
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with the projected propagator D given by 

= D,,^[-A + 2Xa]-'. (5.19) 

The exphcit dependence on a special point xq in Eqs. (5.18), (5.19) resuhs in a nontrivial 
spatial dependence for the solution fields fix, ax- For x far from the fixed spin at xo, 
we expect the field ax to approach the constant value a = —-^^ corresponding to the 
negative dynamically generated squared-mass in the translationally invariant gauge. As 
fixQ is pinned at 1, and fix > 1 va. general, ArZajp > and (from (5.18)) olxq > 0. So the 
saddle-point solution in this case involves a spatially dependent dynamical mass, and 
translational invariance is obviously lost in the propagators D, D above. Nevertheless, if 
we choose periodic boundary conditions at the edge of the lattice, the fixing of a single 
spin still only amounts to a S0(1,A^) gauge fixing, and S0(1,A^) invariant two-point 
correlators such as 

{rix ■ ny)L,N/\,2 ^ rixfiy - XDxy , (5.20) 

must still be translationally invariant, and indeed equal to those found in the transla- 
tionally invariant gauge, namely (5.8). The equations (5.18), (5.19) cannot be solved 
analytically, but they are numerically solvable on a given lattice by iteration: one takes 
a reasonable approximate starting Ansatz for fixi ax and then solves (5.18) for n and a 
in alternation until convergence is reached. Single precision convergence (6-7 digits) is 
typically reached with less than 100 iterations. 

A cross-section of the fi field on a 20 x 20 lattice (with the fixed spin at the center point 
(10, 10)) is displayed in Fig. 10. 
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Figure 10: Cross-section of fla; for L = A = 20. 

The figure exhibits the qualitative features discussed above. The a field solution in this 
case (which corresponds to the parameters of Table 1) has a positive spike at the fixed 
spin with ax^ = +0.5, and ax tending rapidly (within 3 or 4 lattice spacings from the 
fixed spin in all directions) to the negative constant value (recall (5.6)) found in 

the translationally invariant gauge. Note that the typical values of the non-invariant 
quantity fl = {n^)L,N/x,i are very different in the two gauges i = 1, 2. Nevertheless, the 
invariant 2-point functions computed from (5.8) and (5.20) agree, as shown in Fig. 11. 
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Figure 11: Comparison of large N invariant 2-point correlator in fixed spin and translationally 
invariant gauges, L = A = 20. 

6. Conclusions 



We have analyzed nonlinear sigma models with noncompact target space and symme- 
try group S0(1,A'^) in dimensions D > 2 combining analytic and numerical methods. 
The lattice formulation was used with the dynamics defined both in terms of a trans- 
fer operator and a functional integral; in the latter case a gauge fixing was essential. 
Perhaps the most remarkable feature emerging from the analysis is the intricate vac- 
uum structure as witnessed by Theorem 2. Analyzing the system on a finite spatial 
lattice via the transfer matrix - where one would usually expect a unique ground state 
- we could identify a nontrivial ground state orbit, i.e. infinitely many nonnormalizable 
ground states transforming irreducibly under S0(1,A'"). In the thermodynamic limit 
spontaneous symmetry breaking was found to occur in all dimensions D > 1 (where the 
case D = 1 was already treated in [18]). For dimensions less than three this highlights 
that the Mermin- Wagner theorem does not hold for these systems. To (numerically) 
see this spontaneous breakdown on the level of correlation functions, the introduction 
of a suitable new order parameter ('Tanh') was instrumental. The mathematical reason 
for these unusual features was understood to be the non- amenability of the symmetry 
group. 

Since in two dimensions the symmetry breaking is surprising we examined this case 
in more detail. Since the gauge fixing by necessity breaks the symmetry explicitly, 
it was important to study quantitatively the effect of this exphcit breaking via Ward 
identities. Our numerical simulations show clearly that this explicit violation disappears 
in the thermodynamic limit, whereas the symmetry breaking shown by the 'Tanh' order 
parameter remains. The new order parameter thus provides for noncompact models a 
numerically effective way to probe for 'spontaneous' symmetry breaking in finite volume. 

In addition we performed a large N saddle point analysis in the two-dimensional models. 
The qualitative features we deduced for the model at finite N were confirmed explicitly 
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in the solution of the N = oo model. 



A variety of open questions remain: what is the significance of the spontaneous symme- 
try breaking for the localization phenomena the two-dimensional sigma-models provide 
an effective description of? Further are these systems integrable and amenable to a 
bootstrap construction, e.g. based on an i?-matrix with the symmetry of the principal 
unitary series representations like in [54]? For numerical simulations a (hybrid-) cluster 
algorithm would be desirable; in particular to probe whether or not the peculiar long 
range order (sensitivity to boundary conditions) found in the 1-dimensional model ex- 
tends to the field theories. An important question is of course whether there exists a 
nontrivial continuum limit. Conventional wisdom would say no, because the models are 
(perturbatively) asymptotically free in the infrared but not in the ultraviolet ([55, 10]). 
To gain some feeling whether this is true beyond perturbation theory, it might be worth- 
while to study simplified hierarchical versions of the Renormalization Group. Finally 
the vacuum structure in the reconstructed Hilbert space should investigated further, as 
well as its relevance for dimensionally reduced gravity theories, where a continuum limit 
is expected to exist. 
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Appendix A: Spectral decompositions and heat kernel on Mn 



Let — A''^^ be minus the Laplace-Beltrami operator on the hyperboloid H^v, N > 2. 
Recall that its spectrum is absolutely continuous and is given by the interval ^(A^— 1)^ + 
uj'^, uj > 0. There are several complete orthogonal systems of improper eigcnfunctions. 
From a group theoretical viewpoint the most convenient system arc the 'principal plane 
waves' E^^ki'n) (see [56, 25, 27] and the references therein) labeled by a; > and a 
'momentum' vector k E S^"^. Parameterizing n = — 1 s), they read 

KAn) = - Ve^s- k]-'^M-^^ . (A.l) 

The completeness and orthogonality relations take the form 



/ 



[ dujii^ioj) f dS(k)E^An)*E^An') = 5(n,n') , (A.2) 

Jo Js^-^ 

where S(n,n') and 6{k,k') are the normalized delta distributions with respect to the 
invariant measures dfl{n) and dS{k) on Hjv and S^~^, respectively The spectral weight 
is 



^^^^^ ^ J2^ 



r{iuj) 



(A.3) 



The main virtue of these functions is their simple transformation law under SO(l,iV). 
For A e S0\N) one has trivially E^AA''^ri) — E^^Ak{n). To describe the action of the 
boosts let a e S^~^ and decompose n — ^/^^ — 1 s into its components parallel n\\ and 
orthogonal nj_ to a. A boost in the direction a will then leave n±_ invariant. Denoting 
the boost parameter by 6* e M the corresponding element A = A{9, a) = A{—9, a)~^ acts 

by 

^ \ n J \ ch9 n\\ — sh9 ^a + n± J ^ ^ 

Using this in (A.l) one verifies 

Eu,,k{A-^ri) = [ch^ + a • ^sh^]-i(^-^)-^-£;,,,^(fe)(n) , (A.5) 

where fj^{k) G S^~^ is a rotated momentum vector whose components rA{k)\\ parallel 
and orthogonal fA{k)± to a are given by 

^ , r^, a ■ k ch9 + sh.9 ^ 

fA[k)\<^ = a, 

ch6' + k ■ ash^ 

ch6l + A; • ash6l 
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The transformation law (A. 5) characterizes the (even parity) principal unitary series 
Cjv(<^), a; > of S0(1, N), where CAr((^) and its complex conjugate are unitary equiva- 
lent (see e.g. [25], Vol.2, Sections 9.2.1 and 9.2.7). The orthogonality and completeness 
relations (A. 2) amount to the decomposition (3.4) of the quasi-regular representation p 
on L^(E[iv). Further one verifies 



dS{k) = [che + d ■ ksh9f'^dS{rAik)) . (A.7) 

This implies that dS{k) integrals over products of the form E^^kif^)* Et^,k{i^') are invariant 
under the S0(1,A^) action (A. 5). In particular one can define spectral projectors Pi 
commuting with p in terms of their kernels Pi{n ■ n'), Id 

Pi(n-n'):^ I dujpN{oj) I dS{k) E^^k{nyE^^k{n') 
Ji Js^-^ 

JdQ{n')Pi{n ■ n')Pj{n' ■ n") = Pinj{n ■ n") . (A.8) 

Combined with the completeness relation in (A. 2) this shows that the spectra of — A'^'^ 
and of T in (3.5) are absolutely continuous. 

A complete orthogonal set of real cigcnfunctions of — Ahj^ is obtained by taking the 
dS{k) average of the product of E^^^kin) with some spherical harmonics on the /c-sphere. 
This amounts to a decomposition in terms of SO^(A^) irreps where the 'radial' parts of 
the resulting eigenf unctions are given by Legendre functions. Using the normalization 
and the integral representation from ([57] p. 1000) one has in particular 



j^^ jsik) E^M = {2^me - 1)^^^-^) r]_-ZJO . (A.9) 

As a check on the normalizations one can take the ^ — > 1"^ limit in (A.9). The limit on 
the rhs is regular and gives 27r^/^/r(A^/2), which equals the area of S^~^ as required 
by the limit of the Ihs. Denoting the set of real scalar spherical harmonics by Yi^m{k), 
/ e No, m = 0, . . . , d{l) - 1, with d{l) = {21 + N - 2){l + N - 3)\/{l\{N - 2)!) we set 



H^,i,m{n) ■■= J dS{k) Yi^^{k)E^4n) (A.lOa) 

= fc,(^)nm(^i(e^-i)^(^-"^)^iv'^S:\o, with (A.iob) 

//-I nV2 
koicv) = (27r)^/2 , kiicv) = {2nf/' K + + j)"] j ' ^ ^ 1 ■ 

The expression (A.lOb) is manifestly real, the equivalence to (A.lOa) can be seen as 
follows: from (A. 2), (2.3), and the orthogonality and completeness of the spherical 
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harmonics one readily verifies that both (A. 10a) and (A. 10b) satisfy 

[dn{n)H^^l^rn{n)*H^',l',m'{n) = ]^S{UJ - Uj')Si^i>Sm,m' , 

J fiNyi^) 

poo 

/ duj iiN{uj)^H^,i,min)*H^,i,min) = 5{n,n') . (A.ll) 

•^0 l,m 

Further both (A. 10a) and (A. 10b) transform irreducibly with respect to the real d{l) 
dimensional matrix representation of SO^(A'") carried by the spherical harmonics. Hence 
they must coincide. A drawback of the functions (A. 10) is that the k integration spoils 
the simple transformation law (A. 5) under S0(1, A^). The transformation law can now 
be inferred from the addition theorem 

E H^,Un)H^,Un') = i2n f/' [(n-n')^ - l]^^^-^) r'S^^JU^-n') ■ (A.12) 



For example for n' — A-n} this describes the transformation of the SO^(A^) singlet 
H^,ofi{n) under A e S0(1, A^). 

Having laid out the relevant representation theory let us consider the spectral decom- 
position of the various operators under consideration. For the kernel (3.5) of T we use 
an ansatz of the form 

/•oo I* 

t^(n-n';l)=/ dw/iNii^) ><0,n{i^) dS{k) E^^kin)* E^^kin') , (A.13) 
Jo 

chosen such that TE^^^ = ^|3,N{!^)E^J^k■ To determine the eigenvalues we set n' = and 
integrate over k. Using (2.3), (A. 9), and the integral ([57], p. 804) one finds 

A,,..(^) = {2.r'e^D-,:, dUe - l)^^"-^)e-/^«Pi7/^5jO = , (A.14) 

as asserted in (3.7). The spectral representation of the iterated kernel tf}{n-n'; x), x e N, 
equals (A.13) just with A/3^jv(a;) replaced by [A^^jv(<^)]^. 

In view of (3.9), (3.8) this directly yields an integral representation for the heat kernel 
which (after rescaling rg^ jl —>■ r) reads: 

exp(TA«^)(n,n')= / dw 11^(0;) e-^^^'^^'+'^'U dS{k) E^,k{nYE^,u{n') . (A.15) 

Let us briefiy recap the main properties of the heat kernel on Hjv, see e.g. [58] 
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(i) exp(rA'^'^)(n, n') is symmetric in n, n' and is a bi-solution of the heat equation 

(ii) for each n' e H^, dQ{n) exp(TA''^'^)(n, n') is a probabihty measure which converges 
to the Dirac measure S{n, n') as r — >^ 0"*". 

(iii) it is invariant exp{TA^^){An, An') = exp(TA'^^)(n, n'), A e S0{1,N), and hence 
a function of r = arccosh(n ■ n') only, for which we write h^{r). 

(iv) hrij-) is smooth and strictly positive for all r > and r > 0; in particular the 
coincidence limit r — > 0+ is finite. 

Most of these properties are readily verified from the spectral representation (A. 15). 
Properties (i) and (iii) are manifest. The limit lim^^o exp(rA'^^)(n, n') = S{n, n') follows 
from (A. 2). The fact that / rfri(n') exp(rA'^^)(n, n') = 1 for all n G Hat and r > 0, is 
a consequence of (3.9) and (3.6). This gives (ii). The finiteness of the coincidence limit 
is clear from (A. 9) and the remark after it. However the positivity in (iv) is masked by 
the oscillating nature of the Legendre functions. It can be shown, for example, from the 
alternative expressions (A.16), (A. 17) below, where also the smoothness is manifest. 

Using the behavior of the Legendre T'^ functions under a sign flip of ^ one can rewrite 
(A. 15) in a form where a simplified spectral weight appears which equals that of the 
= 1 case for all odd A^ and that of A^ = 2 case for all even A^. With some further 
processing one can show [26] the equivalence to the usual expressions for the heat kernel. 
For A^ = 2 see e.g. [59] Vol. 1, Eqs (3.32), (3.33). For AT > 2 see [26]. The final result 
we quote from [58] 




(A.16) 



K{r) = (27r)-^r-V2e-^^ j 

J r 




) 



2 



,2 



AT even. (A. 17) 



e 4t ^ 



From here one readily verifies the positivity property in (iv) . 
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Appendix B: Finite volume corrections to Ward identities 



Here we derive for Z) = 2 the finite volume corrections to the Ward identity (2.22), 
(2.23). We use the translation invariant gauge fixing 1 of Section 2.1 where the S0(1, N) 
invariance is violated by delta function constraint in Eq. (2.10). One may expect such 
corrections to vanish in the thermodynamic limit. In fact, it is possible to calculate the 
explicit form of these corrections and thereby study directly their volume dependence. 
The most convenient approach starts with a derivation of the exact Ward identity in 
a nonsingular gauge analogous to the A (or ^)-gauges of quantized nonabehan gauge 
theories, and then recovers the delta- function gauge by taking the limit A — > oo. Thus, 
we begin with the functional integral 

Zx - /"j]dn,exp|s'o[n] + ^(^n,)'-^lnn° + 21n^nS|, (B.l) 

X ^ XX X ^ 

with 5*0 [n] regarded as a functional of the 'spatial' components of the n-field only, 
eliminating via n° = Vl + n2. Note that, as in the case of gauge field theory, the 
form of the Faddeev-Popov term is identical in the delta-function and the smooth A 
gauges. 

We shall indicate the procedure for the case of the rotation Ward identity in (2.23) only 
- the derivation of the correction terms for the boost Ward identity is analogous but 
more tedious. We shall comment on it at the end of this Appendix. For the rotation 
Ward identity it is enough to consider a rotation in the n^, r? plane and we may wlog 
consider the S0(l,2) model throughout. Performing a local rotation with angle on 
the n^; field then gives 

"^x ■ "^x+ia. ^ ~ I'^xl I'^x+aI cos {Qx ~ dx+pt + Q^a; ~ 0(x+(i) (B.2a) 

(Xl^^) — ^ (^'^\n^\cos{0^ + a^)^ + (^^\n^\sin{e^ + a^)^ (B.2b) 



4 



:= |fij;| cos^a;, nl :— \nx\sm9x . (B.2c) 



Introducing these transformations into (B.l) and expanding to second order in the ax, 
one generates three sorts of terms, which we shall refer to henceforth as the "A" , "B" 
and "C" type contributions to the rotation Ward identity. The A-terms arise from the 
variation of the jSrix-rix+fi term in the action of (B.l), the B-terms from the A gauge-fixing 
term, and the C terms are those quadratic terms arising as cross terms of the (linear) 
variation in the pure action and gauge-fixing term. Note that the Faddeev-Popov and 
nonlinear field measure terms play no role as they are invariant under rotations. The 
invariance of the functional integral under the change of variables (B.2) then implies 
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A + B + C = where, after some calculation, we find: 

^ = ^ ^ {Ai^a^){A^ay){Ja;i^Jy„) - ^ J^(A^Q;3;)^(ri^ • , (B.3a) 



xj/ xyzw 

C^-pX^{Af,a^)ay{J^i^nlnl) , (B.3c) 



where as in Section 2 



At this point it is convenient to go over to momentum space by introducing discrete 
Fourier transforms appropriate for the lattice in question: thus, = y^p^^^'^^c^ip) 
etc. One then finds for the A type contributions 

^ = E E(2 - 2 cosp,)[2pE' - Mp)] a{p)a{-p) , (B.4) 

with E"" = {n'^n'^^i^) , a ~ 1,2, equal constants by translation and rotation invariance. 
Similarly, the B type terms may be rewritten 

p 

where we have defined two- and four-point functions D and F resp. as 

(n>5) := 6^'Dxy, Dxy = ^Y.^^^'^^''^D{p), 

p 

{nlnlh\0)h\0)) := 1 J] e^^'-(-^)F(p) , 

p 

where the tilde notation indicates Fourier transform on the n fields n"'{p) := e^^'^n^. 
Finally, introducing the three-point function S as follows (A^ denotes a left lattice 
derivative) 

E^y := {{A;jl)nln\0)) = ^J2''^<^~'^E{p) , 

p 

we find that the C type term amounts to 

p 
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To summarize, the exact rotation Ward identity in a A-gauge takes the form 

^(2 -2cosp,)[2(3E' - Mp)] = /3A[2S(p) - D{p) - D{0)] + (3'X'r{p) . (B.7) 

The terms arising from the gauge-fixing part of the action are isolated on the right- 
hand-side of (B.7): the left-hand-side corresponds precisely to the naive Ward-identity 
(2.30). To obtain the form of the Ward identity appropriate for the delta-function 
gauge used for the simulations, we must next examine the limit of the right-hand-side 
when A oo. In order to do this, we note that the action in (B.l) can be written 
S^So + f[n\Of + n^Of]. Furthermore, we have the distributional limit 

e-^-^ S{x) + ^S"{x) + ^,S""{x) + ..., as A ^ oo . (B.8) 

We shall momentarily use the notation {O)o '■= J dn6{J2x^x) Oe~^° to denote the 
expectation of a Green's function G{n) in the delta-gauge, whereas ( • ) will denote the 
A-gauge expectation, as previously. The B term (B.3b) in A-gauge can be rewritten 

X xy 

Using (B.8), it is easy to obtain for the infinite A limit of the first term on the right- 
hand-side of (B.9) 



Again, using (B.8), one easily verifies, for any O independent of n^, in the infinite A 
limit: 

The first derivative term in (B.ll) can be written in terms of a new field ^x, ^ follows: 



dh^{0) V 



with (B.12) 



/ ' 



while the second derivative term is easily seen to be suppressed by a factor of volume 
V relative to the first, and will be ignored henceforth. In momentum space, the B type 
terms are then found to yield 

F E { - ^ + in(^'(p)nH-p)i(Of)} . (B.13) 

P 
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The large A limit of the C type cross term can likewise be evaluated 

/?Aj]a.a,(A;4>Jni(0)) J^(2 - 2 cosp^)a(p)a(-p)L'(p) 

xy P IJ- 

P 

with the modified three-point function 

E'ip) 5]e^^-(--^)(A;j°,nJe(0)) ■ (B.14) 

X 

Combining these results, we find that the rotational Ward identity, given by (B.7) in 
the smooth A gauges, becomes in delta-function gauge 

5^(2 -2cosp^)[2pE' - Mp)] = (B.15) 

-V~V - 2cosp,)D(p) - ^E'ip) + ^,{fi\p)A-p)m') . 

For p fixed, the first two terms arc manifestly of order y for large V. The last two 
terms involve 3- and 4-point functions respectively, with zero-momentum insertions, for 
which the volume dependence is not a-priori clear. However, they may be computed 
easily from the numerical simulations of Section 4. We find that the best fits to the 
volume dependence of the last two terms suggest a behavior ~ In These fits were 
performed using the results of measurements on 32^, 64^ and 128^ lattices. 

Finally, we show in Fig. 12 the right-hand-side of (B.15), divided by the trivial kinematic 
factor ^^{2 — 2 cosp^), for /3 =10 and L =32,64, and 128. In all cases, they represent a 
small numerical correction to the left-hand-side, as expected from the agreement found 
in Section 4. 

Finite volume corrections to the boost Ward identity in (2.23) can be computed in a 
manner precisely analogous to the procedure leading to (B.15). Apart from a trivial y 
term, there are in this case nine structures appearing on the right-hand-side of the Ward 
identity. As the formulas are somewhat lengthy we refrain from spelling them out here. 
However we have studied the finite volume dependence of these terms on 32^, 64^ and 
128^ lattices, and again find that the dominant asymptotic behavior is In V^, as for the 
rotation Ward identity. 
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Figure 12: Volume dependence of correction terms to rotation Ward identity. 
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